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1—5 ; ABSTRACT 

' We present a rapid binary evolution algorithm that enables modelling of even 

the most complex binary systems. In addition to all aspects of single star evolution, 

T-H ■ features such as mass transfer, mass accretion, common-envelope evolution, collisions, 

^ ' supernova kicks and angular momentum loss mechanisms are included. In particular, 

, circularization and synchronization of the orbit by tidal interactions are calculated for 

(N . convective, radiative and degenerate damping mechanisms. We use this algorithm to 

I study the formation and evolution of various binary systems. We also investigate the 

. effect that tidal friction has on the outcome of binary evolution. Using the rapid binary 

' code, we generate a series of large binary populations and evaluate the formation rate 

I of interesting individual species and events. By comparing the results for populations 

■~>^ I with and without tidal friction we quantify the hitherto ignored systematic effect 

' of tides and show that modelling of tidal evolution in binary systems is necessary in 

order to draw accurate conclusions from population synthesis work. Tidal synchronism 

Q ■ is important but because orbits generally circularize before Roche-lobe overflow the 

^ I outcome of the interactions of systems with the same semi-latus rectum is almost 

C/3 , independent of eccentricity. It is not necessary to include a distribution of eccentricities 

' in population synthesis of interacting binaries, however, the initial separations should 

^ I be distributed according to the observed distribution of scmi-latera recta rather than 

• ^ . periods or semi-major axes. 

• Key words: methods: analytic - methods: statistical - stars: evolution - binaries: 

' general - stars: cataclysmic variables ~ galaxies: stellar content 



1 INTRODUCTION 

The evolution of binary stars does not differ from that of 
single stars unless they get in each other's way. If the binary 
orbit is wide enough the individual stars are not affected by 
the presence of a companion so that standard stellar evolu- 
tion theory is all that is required to describe their evolution. 
However if the stars become close they can interact with 
consequences for the evolution and appearance of the stars 
as well as the nature of the orbit. 

The effective gravitational potential in a frame rotating 
with a circular binary system forms equipotential surfaces 
called Roche surfaces. A sphere of the volume enclosed by 
the critical Roche surface defines the Roche-lobe radius of 
each star. If either star fills its Roche-lobe then gas fiows 
from the outer layers of that star through the inner La- 



grangian point that connects the two Roche-lobes. Some or 
all of this gas may be captured by the companion star so 
that mass transfer occurs and, as a result, the subsequent 
evolution of both stars takes a different course from that of 
isolated stars. When the Roche-lobe filling star is a giant, 
with a convective envelope, or is significantly more massive 
than its companion then, as described by Paczyriski (1976), 
the transferred mass may not be captured by the companion 
but instead accumulates in a common envelope surrounding 
both stars. The outcome of common-envelope evolution is 
still not fully understood but possible scenarios include loss 
of the envelope as the two cores spiral-in to form a closer 
binary or coalescence of the two stars. 

Even in a detached system it is still possible for the 
stars to interact tidally. Tides can synchronize the spin of 
the stars with the orbit and circularize an eccentric orbit as 
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the binary tends towards an equilibrium state of minimum 
energy. Further, if one star is losing mass in a stellar wind 
the companion may accrete some of the material with con- 
sequences for the orbit. For a discussion and review of the 
processes involved in close binary evolution and the various 
kinds of binaries or exotic stars that can result see Pringle 
& Wade (1985) and Wijers, Davies & Tout (1996). 

The effects of close binary evolution are observed in 
many systems, such as cataclysmic variables. X-ray binaries 
and Algols, and in the presence of stars such as blue strag- 
glers which cannot be explained by single star evolution. 
While many of the processes involved are not understood 
in detail we do have a qualitative picture of how binaries 
evolve and can hope to construct a model that correctly 
follows them through the various phases of evolution. Ini- 
tial conditions are the mass and composition of the stars, 
the period (or separation) and eccentricity of the orbit. In 
order to conduct statistical studies of complete binary popu- 
lations, i.e. population synthesis, such a model must be able 
to produce any type of binary that is observed in enough 
detail but at the same time be computationally efficient. By 
comparing results from the model with observed populations 
we can enhance our understanding of both binary evolution 
and the initial distributions (e.g. Eggleton, Fitchett & Tout 
1989; Tutukov & Yungelson 1996; Terman, Taam & Savage 
1998; Nelson & Eggleton 2001). 

Models for binary evolution have been presented in the 
past (e.g. Whyte & Eggleton 1985; Pols & Marinus 1994; 
Portegies Zwart & Verbunt 1996). The model we present 
here supersedes the work of Tout et al. (1997) primarily by 
including eccentric orbits and stellar spins, which are sub- 
ject to tidal circularization and synchronization. Amongst 
other improvements the possibility of mass accretion from a 
wind is included. Our model incorporates the detailed single 
star evolution (SSE) formulae of Hurley, Pols & Tout (2000, 
hereinafter Papl) which allow for a wider range of stellar 
types than the description of stellar evolution used by Tout 
et al. (1997). This requires an updating of the treatment 
of processes such as Roche-lobe overflow, common-envelope 
evolution and coalescence by collision. 

Throughout this paper we refer to one star as the pri- 
mary, mass Ail, stellar type fci etc., and the other as the 
secondary (or companion), mass AI2 and stellar type ^2- At 
any time the primary is the star filling, or closest to filling, 
its Roche-lobe. Numerical values of mass, luminosity and ra- 
dius are in solar units unless indicated otherwise. The stellar 
types correspond to the evolutionary phases desiginated by 
the rapid SSE algorithm of Papl, which are: 



— MS star M ^ 0.7 deeply or fully convective 

1 = MS star M ^ 0.7 

2 — Hertzsprung Gap (HG) 

3 = First Giant Branch (GB) 

4 = Core Hehum Burning (CHeB) 

5 — Early Asymptotic Giant Branch (EAGB) 

6 = Thermally Pulsing AGB (TPAGB) 

7 = Naked Hehum Star MS (HeMS) 

8 — Naked Helium Star Hertzsprung Gap (HeHG) 

9 = Naked Hehum Star Giant Branch (HeGB) 



10 = Helium White Dwarf (HeWD) 

11 = Carbon/Oxygen White Dwarf (COWD) 

12 = Oxygen/Neon White Dwarf (ONeWD) 

13 = Neutron Star (NS) 

14 = Black Hole (BH) 

15 — massless remnant. 

The SSE algorithm provides the stellar luminosity L, 
radius R, core mass Mc, core radius Rc, and spin frequency 



for each of the component stars as they evolve. A 



prescription for mass loss from stellar winds is included in 
the SSE algorithm. The algorithm covers all the evolution 
phases from the zero-age main-sequence (ZAMS) , up to and 
including the remnant stages, and is valid for all masses in 
the range 0.1 to 100 Mq and metallicities from Z — 10~* to 
Z — 0.03. This rapid binary evolution algorithm is a nat- 
ural extension of the SSE algorithm. Many of the formulae 
and much of the terminology contained in Papl are utilised 
in this current paper and therefore the interested reader is 
encouraged to review Papl. 

In Section 2 we describe the binary evolution algorithm 
in detail. Section 3 contains illustrative examples of Algol 
and cataclysmic variable evolution and also compares our 
model with the results of certain binary cases highlighted by 
Tout et al. (1997) and other authors. We present the results 
of population synthesis to examine the effects of tides on 
binary evolution in Section 4 and conclusions are given in 
Section 5. 



2 BINARY EVOLUTION 

Before describing our treatment of Roche-lobe overflow 
(RLOF) in detail we describe interaction in detached sys- 
tems, via wind accretion and tides. We also consider how 
mass and angular momentum loss from an individual com- 
ponent affects the system. 

Two stars, bound through their mutual gravity, move 
in elliptical orbits about their common centre-of-mass. In 
plane polar coordinates, r, the separation of the stars, and 
6, the phase angle, the equations of motion for an elliptical 
orbit are 



all 



1 + e cos ( 



and 



h : 



(1) 



(2) 



where a is the semi-major axis of the ellipse and e the ec- 
centricity. The specific angular momentum of the system, 
h = |h|, is given by 



h = r X V 



(3) 



where both r and v = r lie in the orbital plane. Note that 
the semi-latus rectum 



/ = a (1 - e^) = 
where 

Mb = Mi+ M2 



GMb 



(4) 



(5) 
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is the total mass of the system, is constant if orbital angular 
momentum is conserved. These equations consider the stars 
as point masses interacting by gravity alone. Perturbing ef- 
fects, such as tidal forces, are not taken into account. 



2.1 Wind Accretion 

If the primary star loses mass in a wind at a rate Miw, which 
is determined according to a prescription for single stars such 
as that given in Papl, the secondary can accrete some of the 
material as it orbits through it. Wind accretion is important 
for the evolution of ("-Aurigae and VV Cephei systems (Che- 
Bohnenstengel & Reimers 1986; Hack et al. 1992), symbiotic 
stars (Kenyon 1986; Iben & Tutukov 1996) and massive X- 
ray binaries (Lamers, van den Heuvel & Petterson 1976). 

The mean accretion rate, on to the secondary, can be 
estimated according to a Bondi & Hoyle (1944) mechanism 
to be 

2 



(M2a) = 



GM2 



Qw 



2a- 



Miw 



where 



2 



(6) 

(7) 

(8) 



The wind velocity is difficult to determine accurately from 
observations. We set it proportional to the escape velocity 
from the surface, 

^w = 2/3w^, (9) 
til 

where the value of /3w must depend on the spectral type. 
Observations of luminous stars indicate that /3w decreases 
from about 7 for O stars to about 0.5 for A and F stars 
(Lamers, Snow & Lindholm 1995). Cool super-giant wind 
velocities are observed to be 5 — 35kms~^ (Kucinskas 1999) 
where the lower limit should roughly correspond to the wind 
from the largest stars, about 900 R0 . This suggests /3w — 
1/8. 

Averages over an orbital period are justified because 
the duration of mass loss is typically much greater than P. 
Eq. (^ is strictly only valid under the fast wind assumption 
(fw S> Worb). If the orbital and wind velocities are compa- 
rable in size then the binary motion disturbs the shape of 
the wind and the assumption of spherical symmetry is vio- 
lated. The parameter aw is taken to be 3/2 appropriate for 
Bondi-Hoyle accretion and agrees with the lower limit found 
by Boffin & Jorissen (1988) when modelling wind accretion 
from a supergiant to explain the Barium star Capricorni. 

It is necessary that the secondary does not accrete more 
mass than is lost by the primary, as may follow from eq. (^ 
for highly eccentric orbits, so we enforce the condition 



|A/2a| < 0.81A'/iw1 



(10) 



to ensure this, an d for orbital angular momentum considera- 
tions (see Section 2.2). Use of this condition is rare because 



we expect any system with appreciable mass accreti on t o 
also be experiencing strong tidal evolution (see Section 2.3). 
We note that there is no special reason for the choice of 0.8 
as the upper limit. If the mass-ratio of the binary system 



is close to unity it is possible for both stars to lose mass at 
the same time and the primary may accrete some mass from 
the secondary. In this case any interaction between the wind 
material from the two stars is ignored and eq. (^ is simply 
repeated with the roles of the stars reversed. 

If a star loses mass in a stellar wind then angular mo- 
mentum is lost too. When the companion accretes some of 
the material a fraction of the angular momentum lost from 
the the intrinsic spin of the primary goes into the spin an- 
gular momentum of the companion. Therefore 



2 • 2 • 

Jspinl = -Afiw/ll + -^wAflA/i2 



(11) 



where Afiw < and Mia > (see also Section 7.2 of Papl). 
We assume that the specific angular momentum in the wind 
is transferred with perfect efficiency, i.e. pw = 1, whereas 
Ruffert (1999), who performed a numerical 3-D study of 
wind accretion using a high resolution hydrodynamic code, 
found that the efficiency is roughly between 0-70%, depend- 
ing on model parameters. 

Motivated by observations of RS CVn binaries in which 
the more evolved star becomes the less massive before 
Roche- lobe overflow has begun. Tout & Eggleton (1988) sug- 
gested that mass loss may be tidally enhanced by the pres- 
ence of a moderately close companion. They give a simple 
descriptive formula 



M = Mr <^ H- Bw max 



2' RlJ 



(12) 



for the enhanced mass-loss rate where Mn is the Reimers 
rate (cf. Papl), i?L is the Roche- lobe radius and Bw — 10**. 
It is uncertain whether such enhanced mass loss rates are 
realistic for a wide range of binary systems, so we include 
eq. (|l^) in the model with Bw as a variable parameter, but 
typically Bw = is used. 



2.2 Orbital Changes due to Mass Variations 

The orbital angular momentum can be expressed as 



(13) 



which naively gives 

iorb = [(A/iw + Af2A) aj + (Khw + Mia) al] fiorb (14) 

so that if the primary loses mass in a wind the orbit 
loses an amount AMiafflorh of angular momentum but if 
the secondary accretes some of this mass then an amount 
AM2aif2orb of angular momentum is returned to the orbit. 
This is exactly the case of conservative mass exchange where 
the momentum gained by the secondary must be taken from 
the primary. However if the fast wind approximation is con- 
sidered more carefully (Gair & Snellgrove private communi- 
cation, for details see Hurley 2000) then, averaged over an 
orbit 

1 1 



and 
a 



SMi 
' Mb 



Mb ^ 2M2 



2-e" 



Afa 



A/b 



{SM2) 

1 - e2 



(15) 



(16) 



The change to the orbital angular momentum is then 
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{5J) _ 5M1M2 - {5M2) Ml 
~ ^ Ml Mb 

and 



where 



Km-J ^^''-[wJ Ml 



{SM2) 



d i iorb 



2 
0.2 



fMiy 2 
Km-J " 



2 fAhy 2 

= Im: j " 

so that 



Jorb = (^MlWOl ;^M2Aa2 ] fiorb ■ 



Ml 



(17) 



(18) 



(19) 



(20) 



This breaks down if IAM2I |AMi| which should only 
occur in extreme cases, e.g. when R ~ i?L so that t^orb ~ Vvj. 
The restriction to eq. (^ avoids this problem. If mass loss 
from both stars is considered then 



J orb 



/ • Ml • \ o / ■ Mo ■ 

(^A/iw - J^MiaJ at + (^M2w - -J^Mi 



In an eccentric orbit this model takes into account the fact 
that more matter is accreted at periastron than apastron 
so that mass accretion results in circularization of the or- 
bit. The secondary loses momentum owing to the drag it 
experiences in moving through the primary wind and it is 
assumed that these losses are taken from the system in the 
wind. The model gives e/e < in all cases. BofRn & Joris- 
sen (1988) predicted a modest increase in eccentricity, based 
on Huang (1956), due to neglecting changes to the orbital 
velocity. This was subsequently corrected by Theuns, Boffin 
& Jorissen (1996) although their treatment still allows ec- 
centricity growth in certain cases. Another difference stems 
from the time averages where Huang seems to assume, in- 
correctly, that the time average operator is multiplicative. 
Eq. ( |l5[ ) gives a circularization timescale which, for small 
eccentricity, is 

^ = M^,,^f^y^ (22) 



Mb 

where 52 = M2/M1 is the mass-ratio of the secondary star. 
For a 2.0 Mq star on the giant branch the timescale is 



1.263 X 10' 



(23) 



and we show in Section 2.3 that the equivalent timescale due 
to damping of tides raised on the surface of the 2.0 M© star 
is 



0.43 



yr- 



(24) 



Therefore the changes in eccentricity owing to mass varia- 
tions are very small compared to the orbital changes brought 
about tidal friction. Even on the asymptotic giant branch, 
where the mass loss rates increase significantly, this is still 
true because the deep convective envelope even more rapidly 
damps the tides. 



2.3 Tidal Evolution 

Observations of close binary stars reveal that stellar rotation 
tends to synchronize with the orbital motion and can be 
slower (Levato 1974) or faster (Strassmeier 1996) than for 



single stars of the same type. The stars do not need to be in 
contact for this co-rotation to be achieved so a mechanism 
such as tidal friction must be the cause. Tidal interaction is 
important in changing the orbit of a close detached binary. 
The degree of interaction is critically dependent on the ratio 
of the stellar radius to the separation of the stars (Zahn 1977; 
Hut 1981). The presence of a companion introduces a tidal, 
or differential, force which acts to elongate the star along the 
line between the centres of mass, producing tidal bulges. If 
the rotational period of the star is shorter than the orbital 
period then frictional forces on the surface of the star will 
drag the bulge axis ahead of the line of centres. 

The tidal field can be decomposed into two parts: the 
equilibrium tide which, if all forms of dissipation are ne- 
glected, is described by assuming that the star is always in 
hydrostatic equilibrium and the dynamical tide which arises 
from stellar oscillations. If dissipative processes are at work 
within the star the equilibrium tide will lag or precede the 
^2ijine of centres. The resulting torque transfers angular mo- 
mentum between the stellar spin and the orbit. Energy is 
dissipated in the tides which diminishes the total energy. 
Thus the orbital parameters change and either asymptoti- 
cally approach an equilibrium state or lead to an accelerated 
spiralling in of the two stars (Hut 1980). The equilibrium 
state is characterized by co-rotation and a circular orbit, 
corresponding to a minimum energy for a given total angu- 
lar momentum and alignment of the spin-orbit axes (which 
we do not consider in detail). 

The main difficulty in a treatment of tidal evolution 
arises in identifying the physical processes which are respon- 
sible for the tidal torque, that is, the dissipation mechanisms 
which cause the tides to deviate from an instantaneous equi- 
librium shape and thus to be misaligned with the line of cen- 
tres (see Zahn 1992 and Tassoul & Tassoul 1996 for recent 
reviews of tidal evolution theory). 



2.3.1 The Equilibrium Tide with Convective Damping 

The most efficient form of dissipation which may operate 
on the equilibrium tide is turbulent viscosity in the convec- 
tive regions of a star. All other forms of dissipation produce 
timescales that normally exceed stellar nuclear lifetimes. 
This mechanism provides a satisfactory interpretation of the 
behaviour of stars possessing relatively deep convective en- 
velopes. However, it should be noted that this is only an 
interpretation because no completely satisfactory descrip- 
tion of stellar convection is available. Nevertheless, estimates 
based on the eddy-viscosity treatment of convection, i.e. the 
mixing-length theory, seem to adequately represent the tidal 
timescales for such stars (Zahn 1989). 

To investigate the change of binary parameters due to 
tidal friction Hut (1981) considers a simple model in which 
only equilibrium tides are described, with very small devia- 
tions in position and amplitude with respect to the equipo- 
tential surfaces. The resulting tidal evolution equations for 
the tide raised on a star of mass AI due to the presence of 
its companion with mass m are 



de 
dt 



(25) 
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dt 



= 3 



1) slf^) 



(l-e2) 

2\3/2 



(26) 



where the /„'s are polynomial expressions in given by 
Hut (1981). The apsidal motion constant of the primary 
k takes into account the structure of the star and T is the 
timescale on which significant changes in the orbit take place 
through tidal evolution, i.e. the damping timescale defined 
in terms of the time it takes for the tidal bulge to catch-up 
with the current position of the line of centres. Furthermore 
q = m/M (= 52 as defined in this work) is the mass-ratio 
of the stars and rg is the radius of gyration where rg — 
I /{MB?). This set of equations are valid for any value of 
the eccentricity and further discussion of the weak friction 
model can be found in Alexander (1973) and Kopal (1978). 

From eq. ( |2^ ) the synchronization timescale can be de- 
fined as 



(^spin ^orb ) 



and from eq. ( ^5| ) the circularization timescale is 
1 21 / fc\ , fR\^ 



(27) 



(28) 



when e « and, in the case of Tr.\ 



when Sis 



Because a > R, Tsync < Tcirc and co-rotation is achieved be- 
fore the orbit circularizes. When the orbit is nearly circular, 
e ~ 0, it is stable (e < 0) if f2spin/florb < 18/11, in agree- 
ment with Zahn (1977). If fiapin > fiorb then napin < and 
d > which is as expected, i.e. transfer of angular momen- 
tum from the star to the orbit. 

In this work it is more convenient to use the tidal circu- 
larization timescale in the form given by Rasio et al. (1996), 



./"conv -^'-^ci 



a2(l + g2)(^) 



'^circ '^conv 

SO we use 
k 



M 

2 fcoj'w -/V/t 



21 Tc 

in Hut's equations. Here 



cnv — 1 



M 



: 0.4311 



cnv 



1/3 



3L 



yr 



(29) 
(30) 

(31) 



is the eddy turnover timescale (timescale on which the 
largest convective cells turnover) and -Rcnv is the depth of 
the convective envelope. Eq. (pll) is essentially the same as 
eq. (4) of Rasio et al. (1996) except that we use the radius 
in the middle of the convective zone, i.e. r ~ (_R — |_Renv), 
rather than at the base, as the typical position of a convec- 
tive element. The numerical factor /conv is 



/c. 



(^'(2rconv) ) 



where Ptid is the tidal pumping timescale given by 
1 



1 



1 

-Porh 



(33) 



As noted by Rasio et al. (1996), theoretically and obser- 
vationally, /conv = 1 as long as Tconv ^ Poih- However, if 



Ptid < Tconv then the largest convective cells can no longer 
contribute to viscosity, because the velocity field they are 
damping changes direction before they can transfer momen- 
tum. Only the eddies that turn-over in a time less than the 
pumping timescale contribute so that the average length and 
the velocity of these cells are both smaller by the same fac- 
tor 2rconv/ftid. The factor of two arises because tides come 
around twice in each period. 

We implement tidal evolution during a time-step At by 
calculating e and f2spin from eqs. ( ^ ) and ( |26| ) and ensuring 
that the star is not spun down (or up) past the equilibrium 
spin at which no angular momentum can be transferred. 



= ,f2 



Slo 



n3/2 



(34) 



./5(e2)(l 

The rate of change of the rotational angular momentum is 
then given by 

iapin = [fca {M - A4) i?' + fcJ,Mc7?c'] Siapin , (35) 

where k'2 = 0.1 and k'j — 0.21 (see Section 7.2 of Papl which 
contains a description of how Jspin and Jlgpin are calculated) , 
so that 



J orb 
e/anin 



Jorh Jspin At 



Js, 



111 + JspinAt 

eAt 



and the change in semi-major axis is taken care of by the 
change in orbital angular momentum. A description of this 
type has been used successfully for red giants by Karakas, 
Tout & Lattanzio (2000) to model the eccentricities of Bar- 
ium stars. 

The procedure for tidal circularization and synchroniza- 
tion given above is valid for all stars with convective en- 
velopes. These have stellar types 

k = 1 & M < 1.25 
k e {0,2,3,5,6,8,9} . 

If the companion has an appreciable convective envelope too 
then the roles of the two stars can be reversed and the pro- 
cess repeated. Details of how to obtain the mass of the con- 
vective envelope can be found in Section 7.2 of Papl. The 
radial extent of the convective envelope is calculated in a 
similar fashion. For all giant-like stars (fc G {3,5,6,8,9}) 
Ranv — R — Rc- For main-sequence stars 



Ra 




1.25-M \ 1/^ 
0.9 J 



M < 0.35 

0.35 < M < 1.25 

M > 1.25 



(32) ^ = J_ 

tms 



and then 

-Renv = -Rcnv,0 (1 — t)"'"''* 

where 

t 



(36) 



(37) 



(38) 



and R' is the radius of a MS star with M = 0.35 M0 at r (see 
Papl). For Hertzsprung gap stars the convective envelope 
grows as the stars evolve so that 



Rcnv 

with 



r^/' (i?-i?c) 



(39) 
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t — tms 



(40) 



The main-sequence lifetime, tyis, and the time taken for a 
star to reach the base of the GB, tsGB, are given by eqns. (4) 
and (5) of Papl. 

2.3.2 The Dynamical Tide with Radiative Damping 

For stars with radiative envelopes, another mechanism is 
required to explain the observed synchronization of close 
binaries. A star can experience a range of oscillations that 
arise from, and are driven by, the tidal field: the dynamical 
tide (Zahn 1975). 

Zahn (1977) derived a circularization timescale 



1 _ 21 



2 \ ) 



92 (1 + 92 



.11/6 



E2 



^.21/2 



(41) 



for radiative damping of the dynamical tide. Comparison 
with eq. (Eq) gives 



1.9782 X 10*^ (l + q2)-'/«i?2yr-^ 



(42) 



where E2 is a second-order tidal coefficient which can be 
fitted to values given by Zahn (1975), 



E2 = 1.592 X 10"'' M^-*"* . 



(43) 



We use (k/T)^, along with = k'2, in Hut's equations to 
model tides raised on stars with a radiative envelope. These 
have types 

fc = 1 M > 1.25 
k G {4, 7} . 

The corresponding synchronization timescale is given by 



= 52' 



5/3 



/GMy/^ MR^ 



'^synch 

SO that 



ql (1 + 92)'/' E2 



(44) 



"^circ '^synch 



K ^ 0.075 



1 + 9 2 (R 
92 



(45) 



and once again Tcirc > Tsynch for reaUstic mass-ratios. 



2.3.3 Tides on Degenerate Stars 

Campbell (1984) examined the tidal effects in double- 
degenerate binaries and derived a synchronization timescale 



(46) 



Tsynch 1.3 X lO'^ 

for a tide raised on a WD of mass M. This assumes that 
the WD on which the tides act is initially non-rotating in 
inertial space and Tsynch provides an upper limit if this is not 
true. The synchronization timescale would also be shorter if 
tidally excited non-radial oscillations in the star were con- 
sidered. For WD- WD binaries the orbit should already be 
circular and a circularization timescale is not relevant. How- 
ever a synchronization timescale is applicable because the 
companion may be spun-up by mass transfer. The formu- 
lation also applies to WD-NS binaries which are in most 
cases initially eccentric so that a circularization mechanism 
is required. 



M 


-RzAMS 


<Ms/yr 


(a/ -R)sync 


(a/i?)circ 


0.5 


0.46 


1.29 X 10^^ 


125.46 


32.76 


0.8 


0.73 


2.65 X 10"'^'^ 


109.81 


24.64 


1.0 


0.89 


1.10 X 10"' 


107.61 


20.36 


1.2 


1.14 


5.62 X 10^ 


124.58 


14.71 


1.6 


1.48 


2.25 X 10^ 


6.80 


3.94 


2.0 


1.61 


1.16 X 10^ 


6.75 


3.92 


3.2 


2.05 


3.18 X 10** 


6.69 


3.89 


5.0 


2.64 


1.04 X 10* 


6.68 


3.88 


7.0 


3.20 


4.89 X 10'' 


6.74 


3.91 


10.0 


3.94 


2.43 X 10'' 


6.89 


3.98 



Table 1. The fractional separations at which the synchroniza- 
tion and circularization timescales are equal to 1/4 of the main- 
sequence lifetime for equal-mass binaries starting on the ZAMS. 
Only the contribution from one star has been taken into account. 
Column 2 gives the ZAMS radius of the star. 



Comparing the synchronization timescale with that 
from eq. ( p^ ) for convective damping gives 

k\ ._«,/Lx5/^ 



(I), 



2.564 X 10^ 



yr 



(47) 



which can be used along with Tg = ^3 in Hut's equations to 
give an estimate of the effect of tidal forces on the orbital 
parameters owing to degenerate damping. 

In most cases a WD primary is expected to be effec- 
tively non-rotating because the progenitor star would have 
spun-down considerably on the AGB through mass loss and 
a large radius increase. However, it should be noted that 
if theprimary is initially rotating it may be better to take 
eq. (kffl as representative of the circularization time. As- 
suming similar results to those derived for convective and 
radiative damping, this would mean an underestimate of 
Tcirc. Therefore as eq. (^6|) is an overestimate of Tsynch, and 
in the absence of a better understanding, we can hope that 
one approximation corrects for the other. 

2.3.^ Relative Strengths of the Damping Mechanisms 

Timescales for the synchronization and circularization of bi- 
nary stars with radiative envelopes are generally orders of 
magnitude larger than those characterizing the equilibrium 
tide in convective-envelope stars. As an illustration consider 
a binary consisting of two MS stars, with masses 1.0 M© 
and 2.0 respectively, separated by a distance of 10.0 Rq. 
The circularization timescale owing to tides raised on the 

1.0 Mr star is Tcirc 



2.5 X 10 yr as given by eq. 



,mp- 



for convective damping while that due to radiative dam 
ing of the tides on the 2.0 Mq star is, according to eq. ( 
Tcirc, r — 1.8 X 10^" yr, an order of magnitude greater than 
the nuclear lifetime of the star. Tidal friction on the 1.0 M© 
star dominates changes to the orbital parameters because 
the convective timescale, Tcirc, c, is an order of magnitude 
less than the MS lifetime of the 2.0 M0 star. For compar- 
ison a WD-NS binary with a separation 10 times the WD 
radius would have Tcirc, d ^ lO^** yr. 

Comparing the circularization timescales for a binary 
consisting of a 0.5 Mq MS star and a 0.8 M© WD gives 

Tcirc, c ~ 6O(a/R0)*yr 
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'^circ.d 



10^^ {a/Ref yr 



and thus the degenerate damping is only dominant for WD- 
WD and WD-NS systems in which the separation can be- 
come very small. 

Following Zahn (1977) the limiting separations for syn- 
chronization and circularization can be defined as the initial 
separations {a/R)sync and {a/R)circ for which the charac- 
teristic timescales are equal. As an example consider equal- 
mass MS binaries with both Tsync and rdrc set equal to one- 
quarter of the MS lifetime. The results are given in Table |l| 
The tidal damping is clearly more efficient for stars with 
convective envelopes and only becomes more so if the bi- 
nary includes a giant-like star. 

We note that Zahn (1977) obtains larger values of 
(a/_R)circ for radiative stars than those given in Table ^. This 
seems to arise from using MS lifetimes up to six times greater 
than those found here. However, using his {a/R)ciic with 
Tsync ~ Tcirc givcs larger (a/_R)sync values than he records 
with the end result that, while the values of (a/_R)circ cal- 
culated here are lower than his, the values of {a/R)sync are 
much the same. 



2.4 Gravitational Radiation and Magnetic 
Braking 

Coalescing neutron stars are a source of gravitational waves 
and may produce 7-ray bursts (Hartmann 1996). Such merg- 
ing requires angular momentum loss from the system. This 
can be achieved by gravitational radiation from close bi- 
nary systems driving the system to a mass transfer state, 
possibly followed by coalescence. It is also necessary for the 
formation of AM CVn systems (Gonzalez Perez 1999) and 
adequate for cataclysmic variables (CVs) below the period 
gap (Faulkner 1971). Eggleton (private communication) uses 
the weak-field approximation of general relativity to derive 
the orbital changes due to gravitational radiation from two 
point masses as 



J or 



.315 X 10 



yr 



(1 



n5/2 



-8.315 X 10 



(1 



(48) 
(49) 



These formulae do not take into account any contribu- 
tion from distortions in spherical symmetry of the stars. 
Kuznetsov et al. (1998) show that for a binary consisting 
of two neutron stars this contribution only becomes notice- 
able when the separation decreases to less than seven times 
the size of the neutron star. Neglecting it then causes less 
than a 10% error in the small time remaining to coalescence. 

Gravitational radiation is only efficient for CVs with 
orbital periods less than 3 hours and so cannot explain the 
mass transfer rates of such objects with periods up to 10 
hours (Faulkner 1971; ZangriUi, Tout & Bianchini 1997). 
However, it is possible for orbital angular momentum to be 
lost from the system via magnetic braking of the tidally cou- 
pled primary by its own magnetic wind. Rappaport, Verbunt 
& Joss (1983) used mass-transfer rates deduced from obser- 
vations of CVs, and the Skumanich (1972) braking law, to 
parameterize the rate of angular momentum loss by mag- 
netic braking. We use 



Jspin = Jmb = -5.83x10 



;Mc, 



M 



(ijfi. 



sp.„/ M0R^yr-^(5O) 



where all masses and the radius are in solar units and 57spin 
in units of years, to alter the spin angular momentum of the 
component star and allow tides to alter Jorb- This is effec- 
tive for any star with an appreciable convective envelope. 
We note that the non-parametric formulation for magnetic 
braking given by Stgpien (1995), derived and calibrated from 
observations of spin-down of single stars, should be investi- 
gated in future versions of our binary model. 

The main-sequence primary in a CV has Ri ~ -Rli — 
Ml so it is possible to derive the relation P/hr ~ 
IO.OM1/M0 for the period (Verbunt 1984). This with 
eqs. ( pq ) and (jH^) gives 



Jn 



9.815 X 10- 



(MbMi^ 



Mi 



(51) 



so that magnetic braking dominates even as Mi 0.0. Ob- 
servations of CVs reveal a gap between 2 and 3 hr in the oth- 
erwise smooth period-mass distribution. This period gap can 
be explained if magnetic braking becomes less effective when 
the orbital period falls below 3hr, corresponding roughly to 
a main-sequence primary which has just become fully con- 
vective (Rappaport, Verbunt & Joss 1983; ZangriUi, Tout & 
Bianchini 1997). We therefore do not apply magnetic brak- 
ing when the primary is a fully convective main-sequence 
star. Ml < 0.35. Braking proceeds only by gravitational ra- 
diation for such systems. 



2.5 Supernovas Kicks 

When a star explodes as a supernova and becomes a neu- 
tron star or black hole it receives a velocity kick, due to any 
asymmetry in the explosion (Lyne & Lorimer 1994). This 
may disrupt the binary. Evidence for such a kick includes 
observations of pulsars with velocities in excess of those of 
ordinary stars (Hansen & Phinney 1997), double-NS binaries 
(Fryer & Kalogera 1997) and X-ray binaries with large ec- 
centricities (Kaspi et al. 1996). The state of the binary post- 
supernova depends on the orbital parameters at the moment 
of explosion and t he k ick velocity. A complete derivation is 
given in Appendix 



Al 



If the binary survives the supernova explosion then it 
is quite likely that the orbital parameters, a and e, change 
markedly from those of the initial orbit. In addition, the 
mass lost by the primary star carries away an amount 
AJoib ~ AAfiaff2orb of orbital angular momentum from 
the system. Before the primary core collapses it has spin 
angular momentum 



e/spin — k'^KIcRc ^spin 



(52) 



which we conserve in the neutron star, spinning it up as it 
shrinks. 



2.6 Roche-Lobe Overflow 

Mass transfer occurs in close binary systems following the 
onset of Roche-lobe overflow. This can be triggered either by 
a star expanding to fill its Roche-lobe as a result of stellar 
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evolution or by angular momentum losses causing contrac- 
tion of the orbit. The Roche- lobe radius of a star is fitted 
by Eggleton (1983) with 



a 



0.49^1 



2/3 



2/3 



In 1 



1/3 



(53) 



in terms of the semi-major axis of the orbit and the mass- 
ratio qi = M1/M2 of the primary star, accurate to within 
2% for < gi < 00. Note that eq. ( p^ ) may also be used to 
obtain the Roche-lobe radius of the secondary by using the 
appropriate mass-ratio, i.e. 52. The theory of RLOF is based 
on two stars in a circular orbit in which complete co-rotation 
has been achieved. In most cases this is adequate because 
tidal friction generally acts to remove any eccentricity on a 
timescale shorter than the evolution timescale of the binary. 
However it may be possble for RLOF to occur in an eccen- 
tric orbit if the binary is formed by tidal capture so that the 
primary, the more evolved star, has not spent all its life as 
part of the present system. It is also possible that stars may 
form in a close eccentic binary but it is generally expected 
that some initial circularization occurs as part of the forma- 
tion process. Eccentric RLOF could also be envisaged if a 
star in an eccentric orbit is rapidly expanding, such as on 
the AGB, so that the nuclear and circularization timescales 
are similar. If any of these rare cases arise then, for want 
of a more detailed treatment, we subject them to instant 
synchronization at the onset of RLOF. 

When _Ri > _Rli mass is transferred to the companion 
star with the primary losing an amount of mass AMir and 
the secondary gaining some fraction AM2r) of this mass. 
All normal processes of binary evolution, as described in 
the preceding sections of this paper, are also treated during 
RLOF. This means that mass variations due to stellar winds 
and the related changes to Jorb and Jspin are still accounted 
for as are changes resulting from tidal friction. In addition, 
if a degenerate secondary star is not altered by the tidal 
evolution treatment then its spin is changed according to 



AJ,pin2 = k'iAAhRRl fit 



(54) 



where we assume that the accreted material comes from the 
inner edge of an accretion disk which has 



GM2 



(55) 



Care is taken to conserve the total angular momentum of the 
system. Independent of the transfer the tides should main- 
tain co-rotation of the primary and the orbit during the 
Roche phase. We assume that mass changes due to RLOF 
do not affect the orbital angular momentum if the mass 
transfer is conservative. Any material not accepted by the 
secondary is lost from the system taking with it specific an- 
gular momentum from the primary so that a (A/i + M2) = 
constant. During RLOF _Rli is used as the effective radius 
of the primary when calculating Miw and any associated 
angular momentum changes, as well as in the tidal expres- 
sions. 

The treatment of RLOF presented here is a substan- 
tially revised version of that described by Tout et al. (1997) 
including essential changes required to make it compatible 
with the updated stellar evolution treatment and the ex- 
plicit treatment of angular momentum. Following Tout et 



al. (1997) we describe the stability of mass transfer using 
the radius- mass exponents C, defined by Webbink (1985). 



2.6.1 Dynamical Mass Transfer 

If Cad < Cl the radius of the primary increases faster than 
the Roche-lobe on conservative mass transfer. The mass loss 
rate from the primary is limited only by the sonic expansion 
rate of its envelope as mass is transferred through the in- 
ner Lagrangian point connecting the Roche-lobes of the two 
stars. Stars with deep surface convection zones, and degen- 
erate stars, are unstable to such dynamical timescale mass- 
loss unless the primary is rather less massive than the sec- 
ondary so that Cl is more negative than Cad- Thus dynami- 
cal mass transfer occurs for giants (fci € {3, 5, 6, 8, 9}), low- 
mass MS stars (fci = 0) and WDs (fci G {10, 11, 12}) when 
?i > gcrit. This critical mass-ratio is defined by Cad = Cl 
where Cl ~ 2.13gi — 1.67 (Tout et al. 1997) for conservative 
mass transfer. For normal giants R oc , where x ~ 0.3 
varies with composition (as given by eq. (47) of Papl), so 
that Cad = —2; if Cad ~ Ccq- Howevcr the behaviour of the 
radius deviates from this relation when the mass of the giant 
envelope is small. To detailed stellar models we fit 



Cad ~ Cc' 



-x + 2 



M 



and use 



gcrit 



= 1.67 -X + 2 



/MciV 
V Ml J 



/2.13 



(56) 



(57) 



for the critical mass-ratio above which mass transfer pro- 
ceeds on a dynamical timescale. We note that a modified 
criterion is required in cases of non-conservative mass trans- 
fer and that the assumption of Cad ~ Ccq is not always true 
(Hjellming 1989) . Naked helium giants have Cad ~ so that 
Qcrit ~ 0.784. We treat all cases with qi > gcrit in which the 
primary star is a giant as common-envelope evolution, de- 
scribed in Section 



2.7.1 



For low-mass main-sequence stars 
gcrit = 0.695 and for white dwarfs gcrit ~ 0.628 (see Tout 
et al. 1997 ). We describe treatment of these cases in Sec- 



tions 2.6.4 and 2.6.1: 



Based on models of condensed polytropes (Hjellming & 
Webbink 1987) an alternative condition for dynamical mass 
transfer from a giant primary, 

gcrit = 0.362 + [3{1~ Afd/Mi)]"' , 

valid for Md/Mi ^ 0.2, isgiven by Webbink (1988). This 
relation is similar to eq. (pqj for Md/Mi = but quickly 
diverges as Md/AIi increases; it is a factor of two larger 
at Md/Mi = 0.6. For values of qi intermediate between 
the two conditions Tout et al. (1997) allow mass transfer to 
proceed on a thermal timescale so that common-envelope 
evolution is avoided. In the presence of mass loss, and par- 
ticularly any enhanced mass loss, this would lead to a signifi- 
cant increase in the number of stable mass transfer systems. 
The formation paths of many binary populations, such as 
CVs (which require a common-envelope phase) and symbi- 
otic stars (which need to avoid common-envelope), would 
therefore be affected if we were to follow the same scheme. 
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2. 6. 2 Nuclear Mass Transfer 

If Cl < (Cad,Ceq) then mass transfer is stable until nuclear 
evolution causes further expansion of the star, that is to say 
that the mass transfer is not self-stimulating. The radius of 
the primary remains constrained to that of its Roche-lobe 
and the star remains in thermal equilibrium. Mass transfer 
proceeds at a rate such that R\ « i?Li so that the primary 
just overfills its Roche-lobe. We achieve this by transferring 
mass at a rate that steeply increases with the amount by 
which the Roche-lobe is overfilled, 



In this dynamically limited case, if R\ > 107?li as a re- 
sult of significant orbital contraction, the stars are allowed 
to merge ac cordin g to the appropriate collision prescription 
(see Section 2.7.3). 



MlR = 

where 



F(Mi)[ln(i?i/i?Li)]' M0yr-\ 



(58) 



F {Ml) = 3x IQ-'^ [min (Mi, 5.0)]^ 



(59) 

is chosen by experiment to ensure that the mass transfer 
is steady. If the primary is a degenerate star (fci > 10) 
then it is necessary to increase eq. ( [59[ ) by the factor 
10^/ max (_Ri, 10~*) . In general, nuclear mass transfer oc- 
curs on a much longer timescale than either dynamical or 
thermal mass transfer so that systems such as Algols and 
CVs can be observed in this state. 



2.6.3 Thermal Mass Transfer 

If Ccq < Cl < Cad mass transfer is unstable on a thermal 
timescale. The primary does not remain in thermal equilib- 
rium as it loses mass: it contracts and remains just filling its 
Roche-lobe. This is an awkward case to treat because the 
thermal timescale is long compared with the orbital period 
but short compared with the timescale for nuclear evolu- 
tion. First we calculate the mass transfer rate at the full 
equilibrium rate using eq. (^). This gives an upper limit. 
For giants and giant-like stars (fci £ {2,3,4,5,6,8,9}) we 
limit this to the thermal rate which we approximate by 

Ml 



with 
Tkh 

yr 



lO^M 
RL 



M 
(M- 



(60) 



fee {0,1,7, 10,11,12,13, 14} 
k G {2,3,4,5,6,8,9} 



where tkh is the Kelvin-Helmholtz timescale. This still al- 
lows the possibility of _Ri > _Rli by quite a large fraction. 
The most common case of mass transfer at the thermal rate 
has a Hertzsprung gap primary star (fci — 2) because (^ad 
is quite large while ("eq is close to its giant branch value. 
For these stars ("l can still exceed (ad when qi is large. We 
assume (ad ~ 6.85 (Tout et al. 1997) so that if gi exceeds 
Qcrit = 4 common-envelope evolution ensues as a result of 
dynamical mass transfer. We note that this assumption is 
rather approximate and can be improved in future versions 
of the algorithm by calibrating to detailed binary calcula- 
tions. For all non-giant-like stars we allow mass transfer to 
proceed at the rate given by eq. (jHsj) but limited by the 
dynamical rate 

Ml 



''"dyni 



where 



tdyn = 5.05 X 10 




(62) 



(63) 



2.6.4 



Dynamical Mass Transfer from Low-Mass 
Main-Sequence Stars 



Low-mass MS stars (fci — 0) are deeply convective so 
that mass transfer to a companion proceeds dynamically 
if qi > 0.695. If this is the case we assume that the entire 
star overflows its Roche-lobe on a timescale Tm = tdyni and 
only a single star remains. 

If the secondary star is fairly unevolved {k2 € {0, 1, 2}) 
its thermal timescale will be relatively long and it can only 
accrete a fraction 



AM2R = min f 10^^, l) AMm , 

V rKH2 / 



(64) 



of the mass, where AMm = Mi. If the seco ndary is still 
on the MS it must be rejuvenated (see Section 2.6.(: ). When 
the secondary is a giant or a CHeB star (fe G {3,4,5,6}) 
its envelope can easily absorb all of the primary mass on 
a dynamical timescale. If the secondary is a naked He star 
or a WD (fe G {7, 8, 9, 10, 11, 12}) then aU the material is 
accreted dynamically and swells up to form a giant envelope 
around the degenerate star which becomes the core of a new 
giant star. This is effectively a reverse of the process which 
formed the degenerate remnant star in the first place when a 
giant or CHeB star lost its envelope. The new star requires 
an age appropriate to i ts core, the calculation of which is 
described in Section 2.7 along with a detailed description of 
the outcome. 



2. 6. 5 Mass Transfer from Degenerate Objects 

If a WD (fci G {10, 11, 12}) evolves to flU its Roche-lobe then 
the secondary must be a more massive WD or possibly a NS 
or BH (fe G {10, 11, 12, 13, 14} V Mass transfer proceeds at 
a steady rate according to eq. (pq), unless q\ > 0.628 when 
(61) the mass transfer becomes dynamical. In this case the entire 
mass is transferred on a timescale tdyn i , as given by eq. ( |63| ) , 
leading to coalescence of the two stars. 

When dynamical mass transfer leads to the coalescence 
of two HeWDs we assume that the temperature produced is 
enough to ignite the triple-a reaction. If all the material in 
the new core were rapidly converted to CO then the nuclear 
energy released would be greater than the binding energy of 
the core so that no remnant is left (fe2 = 15). Here we follow 
this procedure but it is possible that a naked helium star is 
produced (Webbink 1984). This alternative could easily be 
included. Dynamical mass transfer of helium on to a CO or 
ONeWD causes the accreted material to swell up and form 
an envelope around the CO or ONe core (Iben, Tutukov & 
Yungelson 1996) so that a HeGB star is formed {k2 = 9). 
If an ONeWD accretes CO or ONe material this is simply 
added to the WD (^2 = 12), unless the new mass exceeds the 
Chandrasekhar mass, Mch, in which case electron capture 
on ■^''Mg nuclei leads to an accretion-induced collapse (AIC; 
Nomoto & Kondo 1991; van Paradijs et al. 1997) and the 
formation of a NS {k2 — 13). 

Dynamical mass transfer from a COWD to another 
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COWD quickly leads to the formation of a thick accretion 
disk (as in all the above cases of dynamical transfer) around 
the more massive WD followed by coalescence over a vis- 
cous timescale. When two COWDs coalesce we form a larger 
COWD (^2 — 11) of mass A'h + AAhn- If the mass exceeds 
Mch then it explodes as a possible type la SN (Branch 
1998) leaving no remnant (fe2 = 15). However, whether a 
thermonuclear explosion actually occurs is not clear. The 
temperature produced at the core-disk boundary depends 
on the accretion rate which in turn depends on the viscosity 
of the disk, all of which is uncertain. If the temperature is 
hot enough to ignite carbon and oxygen, whether the WD is 
converted to an ONeWD depends on competition between 
the rate of propagation of the flame inwards, which depends 
on the opacity, and the cooling rate of the WD. Saio & 
Nomoto (1998) used spherically symmetric evolution mod- 
els of accreting WDs to show that carbon burning is ignited 
at the core-disk boundary and that the flame propagates all 
the way to the centre, by heat conduction, converting the 
material to ONeMg without causing an explosion. If this is 
indeed the case then CO-CO WD merger products with a 
total mass greater than Mch would become AIC neutron 
stars rather than type la SNe. 

When a NS or BH (fci G {13, 14}) fills its Roche-lobe 
the secondary must be a NS or BH and we allow the two 
stars to merge and form a single remnant of the combined 
mass. Alternatively coalescing NSs may destroy themselves 
in a 7-ray burst (Hartmann 1996). 



2. 6. 6 The Response of the Secondary 

If the secondary is a MS, HG or CHeB star (fc2 G {0, 1, 2, 4}) 



we transfer an amount of mass given by eq. 
the thermal timescale of the secondary 

Mi 



< TkH2 • 



), limited by 



(65) 



A giant secondary can respond to mass transfer by dynam- 
ically shrinking and no limitation need be applied. When a 
HeMS star accretes mass from a primary of the same type we 
transfer an amount of mass given by eq. (^^. The HeMS star 
must then be rejuvenated. However, if mass is transferred 
to a naked helium star from a primary of any hydrogen-rich 
type, we put all the mass in an envelope around the helium 
core and form a new CHeB or AGB star. Details of how the 



new star is made are given in Section 2.7 



At any time the secondary may respond to mass transfer 
by filling its own Roche-lobe so that a contact binary is 
formed. If t his oc curs we allow the two stars to coalesce 
(see Section 2.1. c). Generally common-envelope evolution 



arises before contact unless the binary formed as a very close 
system. 

Rejuvenation 

We extend the rejuvanation described by Tout et 
al. (1997) for MS stars to include HG and HeMS stars. When 
such a primary transfers mass it must be aged so that the 
fractional age of the star, in terms of its current evolution 
phase, remains unchanged. The process by which this is done 
is the same as that described in Section 7.1 of Papl. If any of 
these stars are accepting mass then they rejuvenate. This we 
do simply by reversing the ageing process for MS stars with 



radiative cores (0.35 < M < 1.25) and for HG secondaries 
(see Section 7.1 of Papl). For MS stars with convective cores 
(M > 1.25), and fully convective stars (M < 0.35), the core 
grows and mixes in unburnt fuel as the star gains mass, so 
that the star appears even younger. The rejuvenation of MS 
stars is detailed in Tout et al. (1997). The core of a HeMS 
star is convective with Mc oc M so these secondaries are 
rejuvenated in the same way. 

Steady Accretion on to Degenerate Objects 

For transfer of mass to degenerate objects we follow 
Tout et al. (1997) but the greater variety of stars now makes 
the procedure more complex. Accretion of H-rich material 
leads to novae if A'/ < 1.03X 10"'^ M0yr~\ and supersoft X- 
ray sources if 1.03x10^'^ < M < 2.71 x lO"'' MQyr"\ in the 
same way and for accretion rates M > 2.71 x 10~^ M© yr~^ 
HeWDs become GB stars while CO and ONeWDs become 



TPAGB stars (see Section 2.7 for details). For nova systems 
we assume that all but a fraction e of the accreted material 
is ejected in the nova explosion, i.e. 



(66) 



where e is an input parameter to the model and may be 
negative. Typically we take e — 0.001. 

How the WD secondary responds to steady accretion 
of He-rich material, either directly (7 < fci < 10) or as 
a consequence of steady hydrogen burning on the surface 
of the WD, depends on its composition. Woosley, Taam & 
Weaver (1986) calculated thermonuclear explosions induced 
by the accretion of helium on to a He WD. Their models 
showed that accretion at 2 x 10~* Mq yr~^ on to a 0.4 M© 
He WD leads to a violent, centrally ignited detonation when 
the mass of the star reaches 0.66 Mq. Nomoto & Sugimoto 
(1977) had previously found detonations at 0.78 Mq. Here a 
HeWD is allowed to accrete material until its mass reaches 
0.7 Mq when it is destroyed (fo = 15). 

A COWD accreting helium can accumulate about 
0.15 Mq of helium-rich material, provided the total mass re- 
mains below Mch, before helium ignites at the base of the ac- 
creted layer. This is an edge-lit detonation (ELD). The det- 
onation front propagates outward through the helium layer 
while an inward propagating pressure wave compresses the 
CO core which ignites off-centre, followed by an outward det- 
onation which destroys the WD (Woosley & Weaver 1994; 
Livne & Arnett 1995). Kawai, Saio & Nomoto (1987), by 
considering steady state spherically accreting models, argue 
that if the accretion rate exceeds 3 x 10~* Mq yr~^ helium 
can burn steadily and non-degenerately. It is not clear that 
such steady burning is possible if accretion is from a disk 
leaving most of the WD photosphere free to radiate. There- 
fore, once a COWD has accumulated 0.15 Mq of accreted 
helium, we set off an ELD which destroys the WD in a type 
la SN. Livne & Arnett (1995) argue that the properties of 
type la SNe are diverse enough to suggest that the progen- 
itors have a range of parameters, such as mass. The more 
luminous events might then correspond to more massive pro- 
genitors. They find that detonations in sub-Chandrasekhar 
mass WDs are a promising explanation for most, if not all, 
type la S Ne. Such SNe are not standard candles (cf. Sec- 



tion 2.6.E): they exhibit considerable diversity in behaviour. 

In the case of steady transfer of He-rich material on to 
an ONcWD, or steady transfer of C-rich material (fci = 
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11, 12) on to a CO or ONeWD we add the material to 
the degenerate core of the secondary. If the new mass of 
a white dwarf secondary exceeds the Chandrasekhar limit, 
Mch, then we annihilate it in a type la SN, unless it is an 
ONeWD {k2 = 12) which undergoes AIC leaving a NS rem- 
nant. When such a SN explodes as a result of steady mass 
transfer the primary survives, having transferred only just 
enough material for the secondary mass to reach Afch- 

The amount of mass that a WD, NS or BH can accrete 
may be limited by the Eddington limit (Cameron & Mock 
1967), 

MEdd = 2.08 X 10"^(1 +X)-^i?2 M0yr-\ 



where X is the hydrogen mass fraction, so that 
AM2R = min (MEdd-rx/, AMm) . 



(67) 



(68) 



There is some uncertainty as to whether the Eddington 
limit should actually be applied because the energy gen- 
erated in excess of the limit might be removed from the 
system in a strong wind or asymmetrically through a disk. 
Super-Eddington accretion rates may be important in the 
formation of low-mass X-ray binaries and millisecond pul- 
sars (Webbink & Kalogera 1997) and in models of X-ray 
emission from quasars due to accretion from a disk on to 
a rapidly rotating BH (Beloborodov 1998). Imposing the 
Eddington limit significantly reduces the formation rate of 
Type la SNe (Livio 2000). For these reasons the Eddington 
limit, eq. (^sj), is only optionally included. A NS may gain 
enough mass to exceed the maximum NS mass of 1.8 Mq 
(see Section 6.2.2 of Papl) in which case it becomes a BH. 



2.7 Common-envelopes, Coalescence and 
Collisions 

As a result of Roche-lobe overflow it is possible for the bi- 
nary components to come into contact and coalesce or for the 
binary to reach a common-envelope (CE) state. The most 
frequent case of common-envelope evolution involves a giant 
transferring mass to a main-sequence star on a dynamical 
timescale. Although the process is difficult to model, and 
therefore uncertain, it is envisaged that the secondary is not 
able to accept the overfiowing material owing to its rela- 
tively long thermal timescale. The giant envelope overfills 
the Roche-lobes of both stars so that the giant core and the 
MS star are contained within a common-envelope. Owing 
to its expansion the envelope rotates slower than the orbit 
of the core and the MS star so that friction causes them to 
spiral together and transfer energy to the envelope. This pro- 
cess may release sufficient energy to drive off the entire enve- 
lope, leaving a close binary containing a WD and a MS star, 
or it may lead to coalescence of the giant core and the MS 
star. Evidence for a mechanism such as common-envelope 
evolution is provided by CVs and close double-degenerate 
binaries whose characteristics can only be understood if a 
significant amount of angular momentum and mass have 
been removed from the precursor system. Planetary neb- 
ulae containing a close binary at their centre are probably 
the result of a recent common-envelope event (Bond, Liller 
& Mannery 1978). 

In contrast to the relatively gentle nature of coalescence 
it is possible that two stars in a close eccentric orbit could 
collide at periastron before either star has a chance to fill 



its Roche-lobe. Furthermore, in a dense environment such 
as the core of a star cluster, two individual stars may be 
involved in a hyperbolic collision. Generally this involves 
formation of a binary via tidal capture and then common- 
envelope evolution or coalescence, owing to the relatively 
low velocity dispersion in star clusters. We assume that any 
collision involving a giant or giant-like star that has a dense 
core and an appreciable envelope leads to common-envelope 
evolution. All other collisions simply result in a merger of 
the two stars. The outcome of a collision depends on the 
impact parameter, which in turn depends on the relative 
velocity and the relative sizes of the stars, as well as their 
structure. Detailed modelling of specific collision cases has 
been undertaken (Rasio & Shapiro 1991; Benz & Hills 1992; 
Bailey & Davies 1999), but in the absence of results covering 
the full parameter space we employ the simple treatment 
given in the following sections. 



2.7.1 Common-Envelope Evolution 

Common-envelope evolution occurs either as a result of 
a collision between a star with a dense core (fci G 
{2,3,4,5,6,8,9}) and any other star, or at the onset of 
Roche-lobe overfiow where mass is transferred from a gi- 
ant (fci £ {2, 3, 5, 6, 8, 9}) on a dynamical timescale and 
qi > (/crit . The primary is therefore a giant or giant-like star 
with core mass Md and core radius Rd, and an envelope 
mass Menvi = Ml — Md. To make the calculations simpler 
we define effective values, marked with a prime ('), for the 
secondary. If the secondary is a MS star (fe G {0, 1, 7}) then 
it has an effective core mass Mc2 = M2, and an effective core 
radius i?^2 = R2, but Mc2 = 0.0. For a giant or giant-like sec- 
ondary (fc2 G {2,3,4,5,6,8,9}) the effective values are the 
actual values, i.e. Mc2 = Mc2 and J?c2 ~ Rc2, and for degen- 
erate secondaries (^2 > 10) where Mc2 = M2 and 7?c2 = R2. 
The effective envelope mass is given by M^nv2 = M2 — M'^2- 
Our treatment of CE evolution follows closely the de- 
scription of Tout et al. (1997) where the outcome depends 
on the initial binding energy of the envelope and the initial 
orbital energy of the two cores. First we calculate the total 
binding energy of the envelope. 



F - ^ 

-fibind.i — — 



Ri 



M2M;, 



R2 



(69) 



with A = 0.5. Note that we use the indices i and / to rep- 
resent initial and final quantities respectively, in relation to 
the CE event. The initial orbital energy of the cores is 



£0 



1 GMdM'^2 



(70) 



where Ui is the semi-major axis at the onset of the CE phase. 
We assume that the cores spiral-in, transferring orbital en- 
ergy to the envelope with an efficency acE, which is neces- 
sarily a free parameter due to uncertainty in its value. It is 
probably not a constant (Regos & Tout 1995) but generally 
QcE ~ 1 is used. The final separation a/, if sufficient en- 
ergy were released to drive away the entire envelope, can be 
calculated via 



-Ebind.i — QCE {Eoih.f — -Borb.i 

where 



(71) 
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^orb,/ — 



2 af 



(72) 



(Ml + Ma) (Menvi + M^n.a) 



(78) 



Eq. (g3|), in conjunction with a/, gives the Roche-lobe of 
each core: -Rli,/ uses qi = Mci/M^a and -Rl2,/ uses ga = 
Mc2/Mci. If neither core fills its Roche-lobe then a binary 
composed of the cores with separation a/ survives and the 
entire envelope escapes the system. The envelope has been 
removed from one, or both, of the stars leaving behind the 
appropriate remnant (see Section 6 of Papl). We assume the 
envelope ejection is isotropic and both cores emerge from 
the CE in co-rotation with the orbit. We note that HG stars 
(fc G {2, 8}) have an evolving density profile that is not yet 
as steep as in giants. Our assumption that HG stars have a 
core mass growing linearly in time (see Papl) mimics this to 
some extent, although the actual remnant mass in a detailed 
computation may still be somewhat different. 

If Roche-lobe overflow would have occurred for either 
core then they could not have spiralled in so far, instead 
coalescing at an earlier time. We assume this takes place 
when = Rl2, or Rd = Rhi depending on which core 
fills its Roche-lobe first, so that the corresponding a — ul 
can be found from eq. (p3) and then the binding energy left 
in the remaining envelope can be calculated according to 



-Ebind,/ — -Ebind.i = QCE 



1 GMclM^ 

2 TTl 



(73) 



This results in formation of a new giant or giant-like star 
with 

fca G {2,3,4,6,8,9} 
M3 = Mf 

Mc3 = M,i + Mc2 , 

unless fca = 7 and fci G {2,3,4,5,6}, in which case Mc3 — 
Mci + M2 (see later). The resultant mass Mf must be de- 
termined. The star has an envelope binding energy 

Sbind,/ = = (74) 

and Rf can be estimated by noting that on a dynamical 
timescale the radius of a giant responds to changes in mass 
according to i? oc M~^ (where x is given by eq. (47) of Papl). 
The radius Ri of the star, if the system were to coalesce 
immediately, satisfies 

G (Ml -f Afa ) (Ml -I- Ma - Mc3 ) 



£^1 



XR^ 



and as CE evolution is over on a dynamical timescale, 



Rl 
Ri 



Ml + Ma 



M 



Therefore 

J^bind.J _ 
-Ebind.i 



( Mf \ 



l + x 



Mf - Mc: 



Ml -I- Ma - Mc: 



(75) 



(76) 



(77) 



which we solve for Mf by a Newton-Raphson iteration. We 
assume the new star rotates with the same period as the 
orbit immediately prior to coalesence. 

A different treatment is described by Iben & Livio 
(1993). Their initial envelope binding energy is not that of 
the giant (or giants) but that of the initial configuration of 
the CE itself, which surrounds the cores with a diameter of 
about 2ai. This amounts to replacing eq. (p9|) with 



(see Yungelson et al. 1994) giving a reduction in the envelope 
binding energy and therefore less likelihood that the process 
ends in coalescence of the cores. A similar effect can be mim- 
icked in our treatment by increasing Uqe to values greater 
than unity. At first sight this may seem unphysical given the 
definition, but as discussed by Iben & Livio (1993), an in- 
crease in acE can be envisaged if additional energy sources 
other than the orbital energy are involved. Processes with 
potential to supply such energy include enhanced nuclear 
burning in shell burning zones of giants, nuclear burning 
on the surface of a degenerate secondary, dynamo genera- 
tion of magnetic fields and recombination of the hydrogen 
and helium ionization zones in giants. Possibly the CE ab- 
sorbs ordinary nuclear energy in the process of swelling up 
to a diameter of 2ai, but this should occur on a thermal 
timescale. Unfortunately, the theoretical determination of 
reliable values for Qce have proven difficult due to a lack of 
understanding of the processes involved, and our ability to 
model them. 



2.7.2 Coalescence of CE Cores 

When two cores spiral into each other the outcome depends 
on their relative density. If they are of similar compactness, 
we assume that they coalesce and mix completely. If one core 
is considerably more compact than the other, that core sinks 
to the centre without mixing while the less dense core mixes 
with the envelope. All cores, including remnant stars, are 
considered to be compact unless the core is a MS star. The 
case of an EAGB star (fc = 5) is particularly interesting 
because it has a carbon core growing within a stationary 
helium core. Detailed models show that the mass inside the 
hydrogen shell is dense enough for the helium core to be 
used as the core for collision purposes. 

In our model the stellar type fcs of the new star pro- 
duced as a result of coalescence is given by the intersection 
of the fcith column and the fcath row of the collision matrix 
in Table The mass of the new star is M3 = Mf and the 
core mass is Mc3 = Md -I- Mca. If, by chance, there is no 
remaining envelope the coalescence product is the remnant 
remaining if a star of type fcs loses its envelope. Note that 
fcs G {13,14} signifies the creation of an unstable Thorne- 
Zytkow object (Thorne & Zytkow 1977) when the merger 
involves a NS or BH. For simplicity our Thorne-Zytkow ob- 
jects rapidly eject their envelopes leaving only the NS or BH 
(Fryer, Benz & Herant 1996) so that M3 = Ma. 

If fei G {2,3,4,5,6} and fca G {0,1} then fes = fci, 
the MS star is simply absorbed into the giant envelope. The 
initial mass and age of the giant are unaffected by the merger 
because the core mass of the star is not changed. Similarly if 
fci G {8, 9} and fca = 7 because the CO core of the giant sinks 
to the centre becoming the core of the new star. If fci G {8, 9} 
and fca G {0, 1} then fcs = 6 because the MS star mixes with 
the envelope of the helium giant to form a hydrogen-helium 
envelope around the CO core of the giant. An initial mass. 
Mo, 3, and age, ts, appropriate for the core mass and stellar 
type of th e new star must be set. We describe our procedure 
in Section 2.7.4 . This must be done in all merger cases where 
the new star is of a different type to either star involved in 
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Table 2. The stellar type of the new star fca produced when two stars merge as a result of a collision (normal type) or common-envelope 
evolution (bold italic type). 



the merger, or if the core mass of the star has changed as a 
resuh of the merger. 

If fci € {2,3} and G {2,3, 10} then fca = 3 because 
the two helium cores merge to form a larger hehum core 
within the existing envelope (or mixed envelopes). If both 
merging cores are degenerate (Mo < MhcF for the giants) 
then enough heat is produced in the merger to ignite the 
triple-a reaction and the nuclear energy released is greater 
than the binding energy of the new core, so ks = 15. 

If fci € {2, 3} and k2 G {6, 11, 12} then fea = 5 because 
the helium core surrounds the CO (or ONe) core within the 
envelope (or mixed envelopes) to form an EAGB star. The 
same occurs if fci = 6 and k2 = 10. 

If fci = 6 and G {6, 8, 9} then the two cores mix to 
form a larger CO core, surrounded by the mixed envelopes, 
so that fca = 6. This is also the result when fci = 6 and 
fc2 G {11, 12} although only one envelope is involved in this 
case. The binding energy of the new envelope must be re- 
duced by the amount of energy liberated in merging the two 
degenerate cores. This removes the envelope in all cases so 
that in actual fact fca G {11, 12}. 

If fci G {8, 9} and fc2 G {8, 9, 11, 12} then fcs = 9 because 
the CO (and possibly neon) cores mix to form a larger core 
surrounded by a helium envelope. When fci G {8, 9} and 
fc2 ~ 10 we assume that the two stars mix to form an evolved 
HeMS star, fca = 7, with Mo 3 = M3 and 



( \, 
VMci-hMa^ 



(79) 



If fci G {2, 3, 4, 5, 6} and fc2 = 7, fca = 4 because the 
giant core material mixes with the HeMS star to form a new 
core that contains some unburnt helium. All the remaining 
cases that arise as a result of common-envelope evolution, 
i.e. fci G {2,3} and fc2 G {4,5,8,9} as weU as fci G {4,5} 
and fc2 G {4,5,6,8,9,10,11,12}, also have fcs = 4 because 
the merged core is a mixture of processed and unprocessed 



helium. The age of the new CHeB star depends on how ad- 
vanced, in terms of central helium burning, the two progen- 
itor stars were. How this is taken into account is the subject 



of Section 2.7.4 



2.7.3 Collision Outcomes 

It now remains to describe the result of all collisions that do 
not proceed via common-envelope evolution. We allow the 
stars to merge without mass-loss, M3 = Mi +AI2 (cf. Bailey 
1999). The stellar type of the merged product is determined 
by the nature of the colliding stars as given by the collision 
matrix (Table ||). 

As in Tout et al. (1997) the collision of two MS stars, 
fci G {0, 1} and fc2 G {0, 1}, leads to ka G {0, 1}, with Mo,3 = 
A/3 and kg G {0, 1}, with Mo,3 = M3 and 



t3 = 0.1 



iMS3 

"Ms" 



Miti _^ M2t2 



iMSl 



tMS2 



(80) 



where ti and t2 are the ages of the colliding stars. When two 
HeMS stars collide we apply the same except that fci = fc2 = 
fc3 = 7 and the tMsi are replaced by inoMsi. If fci G {0,1} 
and k2 — 7 then fca = 4 because the naked helium star 
sinks to the centre forming a CHeB star with Mc3 — M2. A 
HeWD, fci = 10, colliding with a MS star, fc2 G {0, 1}, sinks 
to the centre of the MS star where hydrogen shell burning 
can begin at its surface. The result is a GB star, fca = 3, 
with core mass Mcs — Ah. Similarly if fci G {11,12} and 
fc2 G {0, 1}, double shell burning begins at the WD surface 
so that fca = 6. If fci — 10 and fc2 = 7 the naked helium 
star absorbs the HeWD leaving a rejuvenated HeMS star, 
fca = 7, with Afo,3 = Ma and 



t3 



M2t2 



(81) 



A/3 tHoMS2 

However if fci G {11, 12} and fc2 = 7 then fc3 = 9 because 
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the WD forms the core of an evolved heUum giant with 
Mc3 = Ml. 

When two HeWDs, ki = 10 and k2 = 10, colhde we as- 
sume that the temperature produced is hot enough to ignite 
the triple-a reaction and that the subsequent nucle ar run - 
away destroys the system so that fca = 15 (cf. Section 2.6.5). 
If fci = 10 and k2 £ {11, 12} then fcs = 9, with Mc3 = M2, 
because the helium material swells up to form an envelope 
around the CO (or ONe) core, which in general is smaller 
and denser by virtue of a larger mass. If fci = 11 and ^2 = 11 
then a larger COWD is formed, fca = 11, unless M3 > Mch 
in which case the new star is destroyed in a type la SN and 
fes = 15. Similarly, if fei = 12 and ^2 £ {11, 12} then fca = 12 
unless M3 > Mch causing th e ONe WD to undergo an AIC 
to a NS, fca = 13 (see Section 2.6.5). 

If fci G {13, 14} and k2 G {0, 1, 7} an unstable Thorne- 
Zytkow object results, i.e. fes = k\ and h'h ~ Mi. Finally, if 
fci G {13, 14} and G {10, 11, 12, 13, 14} then ks = ki, un- 
less the new mass of the NS exceeds the maximum allowable 
NS mass of 1.8 M© in which case we collapse it to a BH. 



2. 7.4 The New Star 

When a new giant or core helium burning star is made (fca G 
{3,4,5,6,8,9}) then an age, t^, and an initial mass, Afo,3 
must be assigned to the star appropriate to its type fcs, mass 
M3 and core mass Mcs. 

If ^3 = 3 then our goal is to find Afo.s so as to place the 
star at the base of the GB with its correct core mass. The de- 
cision as to where the star should start its evolution is quite 
arbitrary, with the base of the GB (BGB) chosen because it 
simplifies the process. As discussed in Section 5.2 of Papl the 
helium core of a GB star is degenerate if M < Mhof (where 
MhbF — 2.0 is the maximum initial mass for which He ignites 
degenerately in a helium fiash) and non-degenerate other- 
wise. So the first step is to place MucF into eq. (44) of Papl 
to find the maximum degenerate BGB core mass. If Afc3 is 
greater than this value then Mo, 3 can be fixed by setting 
Mc,BGB = Mc3 and inverting eq. (44) of Papl. Otherwise it 
is necessary to find the luminosity corresponding to Mc3 us- 
ing the GB Mc-L relation given by eq. (37) of Papl. We use 
a Newton- Raphson procedure to find Mo, 3 such that the lu- 
minosity at the BGB (given by eq. (10) of Papl) equals that 
luminosity, taking Mo, 3 = MhoF as an initial guess. Once 
Mo, 3 has been found we set to the corresponding time 
taken for such a star to reach the BGB as given by eq. (4) of 
Papl. We note that stars with initial mass Mo > Mfgb (see 
Section 5 of Papl) do not have a GB which means that a 
GB star cannot be made with a core mass greater than the 
value given by inserting Mfgb into eq. (44) of Papl. Such a 
situation is extremely unlikely, but if it does occur then we 
form an CHeB star. 

We make a CHeB star, kg = 4, whenever the new core is 
composed of a mixture of burnt and unburnt helium and the 
surrounding envelope contains mostly hydrogen. The new 
core is assigned a fractional age which we determine accord- 
ing to the amount of central helium burning that has already 
occurred in the colliding stars. Each of these stars is assigned 



a fraction y where 

r 



t-tr 



t 



^HeMS 

and then 

J/lA/cl + 1/2 Mc2 



k G {0, 1,2,3, 10} 
fcG {6,8,9,11,12} 

k G {4,5} 

k = 7 



yz = 



(82) 



where Mc must be replaced by M for k = 7. The age of the 
new star is then 



ts — tHcl3 + y3iHc3 



(83) 



where tucia and tu,,3 are dependent on Mo, 3. The time of 
central helium ignition, tuei, the lifetime of the CHeB phase, 
tuB, the time taken to reach the second dredge-up phase at 
the start of the TPAGB, tuv, and the MS lifetime of a naked 
helium star, tnoMs, are given by eqns. (43), (57), (70) and 
(79) of Papl. There is a minimum allowed core mass for 
CHeB stars which corresponds to an initial mass of MhbF. 
To find the minimum possible initial mass relevant to Mcs, 
corresponding to j/3 = 1, it is first necessary to find the 
base of the asymptotic giant branch (BAGB) core mass for 
Mo, 3 — MhcF, using eq. (66) of Papl. If Mcs is larger than 
this value we find Mmin by setting Mcs = A'/c,bagb and 
inverting eq. (66) of Papl, otherwise Mmin ~ MhcF . We find 
the maximum possible initial mass corresponding to 2/3 = 
by setting Mc3 — Mc,hcI and solving for Mmax. A bisection 
method can be used to find Mo, 3, where Mmin < Mo, 3 < 
Mmax, and 

Mc3 = Mc,HcI (Mo,3)+t/3 (AfcBAGB (Afo,3) - Mc,HcI (Mo.s)) .(84) 

If hi — 5 we assume that Mc3 is representative of the 
helium core mass, Afc,Ho, of the new EAGB star and find the 
initial mass by setting Mc3 — Mc,bagb and solving eq. (66) 
of Papl for Afo,3. The age of the star is the time taken to 
reach the BAGB, 



(85) 



For a new TPAGB star, ks = 6, Mcs is the CO core mass 
and our aim is to put the star at the start of the TPAGB 
phase so we set Mc3 — Mc,du. To do this it is necessary to 
find the corresponding helium core mass, Mc,hc, just before 
the transition from the EAGB to the TPAGB. Stars with 
0.8 < Mc,He < 2.25 undergo second dredge-up during the 
transition so that 



Afc,DD = 0.44Afc,Hc + 0.448 
while lower mass stars have 



Mc, 



Mc,Ho 



(86) 



(87) 



Stars with Mc,ho > 2.25 on the EAGB would become a NS 
or BH hole before the TPAGB was reached. Thus if Mc3 > 
1.438, given by Mc,hc = 2.25 in eq. (^6[), the star cannot 
exist with ^3 = 6 and we set it up so that it immediately 
becomes a NS or BH by setting Afc3 ~ A/c,sn (see eq. (75) 
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of Papl) in eq. (g). Therefore 

(Mc3 + 0.35) /0.773 M^a > 1.438 
Mc.Hc = < (Mc3 - 0.448) /0.44 0.8 < < 1.438 



Me: 



Mc3 < 0.8 



and Mc,He = Mc,bagb so that Mo, 3 can be found by solving 
eq. (66) of Papl. The age of the star is = tovs by eq. (70) 
of PapL At this stage it is possible that A/0,3 < if Mc3 is 
less than the minimum allowed core mass for the start of the 
TPAGB. This could happen in the rare event of an extremely 
low mass COWD gaining only a very small envelope to make 
it a TPAGB star. If this occurs then we simply add the small 
amount of mass to the WD so that fcs = 11 and Mo, 3 = M3. 

If fcs G {8, 9} we set Mo,3 and so that the new star 
begins life at the end of the naked helium star MS. Mo, 3 
must be found by a bisection method because some of the 
parameters involved depend on the chosen value of Mo, 3. At 
each iteration we find the luminosity at the end of the main- 
sequence, I/THo, by putting r = 1 in eq. (80) of Papl, and 
the HeGB luminosity corresponding to Afc3 , calculated with 
eq. (84) of Papl, for the current guess of Mo,3. We iterate 
until the two values agree to within a difference of 10~*. 
Once Mo, 3 is determined we set the age, ts = tHoMS3 using 
eq. (79) of Papl. 



2.8 The Evolution Algorithm 

The initial state of a binary system is described by the 
masses, A/i and M2, of the component stars, and the pe- 
riod, P, and eccentricity, e, of the orbit. A choice must also 
be made for the metallicity, Z, of the stars. In general the 
system begins with both stars on the ZAMS but it is possi- 
ble to start the binary star evolution (BSE) algorithm with 
the stars in an evolved state. The stars are each assumed 
to have an initial spin frequency on the ZAMS independent 
of the properties of the orbit (given by eq. (108) of Papl), 
however options exist to begin the stars in co-rotation with 
the orbit or with any given value. A number of input pa- 
rameters and options exist for the BSE code, many of which 
have been introduced in the proceeding sections. These are 
summarized in Table ^ along with their ranges and default 
values. 

Each star has its own stellar evolution time-step, Sti, 
(given by eq. (112) of Papl) which includes restrictions that 
prevent the star changing its mass by more than 1%, or its 
radius by more than 10%. For binary evolution it is partic- 
ularly important to keep the time-steps relatively small so 
that the radius of the primary does not increase too much 
over any one time-step. This aids identification of the time 
when the star first fills its Roche-lobe, if indeed it does. We 
also calculate an orbital time-step such that the angular mo- 
mentum of the system changes by 2%, 



= 0.02- 

and take the actual binary evolution time-step when the 
system is in a detached state to be 



Parameter 


Range 


Default 


Section 


/3w 


(0.125 


7.0) 


0.5 


2.1 










. n\ 

— > Z.UJ 


1.5 








Mw 


0.0 


1.0 


1.0 


2.1 






-Bw 


(0.0 - 


-* 10«) 


0.0 


2.1 






€ 


-1.0 


1.0 


0.001 


2.6 


.6 






(0.5 - 




3.0 


11 


.1 




(Tfe km s"'^ 


0.0 - 




190.0 


Al 






Tidal Evolution 


on 


/off 


on 


2.3 






Eddington Limit 


on 


/off 


off 


2.6 


.6 




Stellar Winds 


on 


/off 


on 


7.1 


of Papl 




(0.0 


2.0) 


0.5 


7.1 


of Papl 



A value of zero corresponds to no velocity kick. 

Table 3. Input parameters and options in the BSE algorithm 
set by the user. Column 2 gives the range of values that each 
parameter may assume (suggested range if bracketed) and default 
values are given in Column 3. The relevant section where each 
parameter or option is discussed is given in Column 4 (note that 
some parameters are introduced in Papl). 



the necessary tidal and braking mechanisms as well as evolv- 
ing each star according to the SSE prescription given in 
Papl. If by virtue of the stars growing in size, or the orbital 
separation shrinking, the periastron separation is less than 
the combined radii of the two stars then the stars collide 
and the system is dealt with according to the prescription 



described in Section 2.7 



If the primary fills its Roche-lobe then we interpolate 
within the last time-step until the radius of the star only 
just exceeds its Roche-lobe radius (1 < _Ri/_Rli < 1.002). 
The sys tem enters RLOF and we treat it according to Sec- 



tion 12. 61 First we test for dynamical mass transfer and deal 
with it if necessary, by merging or common-envelope evo- 
lution, otherwise the mass transfer is steady. We set the 
time-step during RLOF to 

where initially 
10"^ 



(90) 



5t- 



(91) 



where 5t is the previous time-step for the detached system. 
After implementing mass transfer, and any mass loss due 
to stellar winds, and adjusting the separation accordingly, 
as well as implementing the tidal and braking mechanisms, 
we check whether the primary still fills its Roche-lobe. If 
-Ri < i?Li then the system leaves RLOF to once more evolve 
as a detached binary. If Ri > i?Li and R2 > Rl2, a contact 
system has formed and we treat this as either a collision or 
common-envelope evolution depending on the nature of the 
stars involved. Otherwise we allow the time-step to grow by 
a factor of two, until Ri/R^i is such that 5Mi/M\ — 0.005, 
i.e. 



5t = min {5ti, 5t2, 5th) 



(89) 



2,ku 



0.005A/1 

\SAh\ 



(92) 



We evolve the system forward by this time, implementing 



and the RLOF process is repeated. 
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Figure 1. Hertzsprung-Russell diagram siiowing the evolution of 
a 2.9 M0 star (* points) and a 0.9 Mq star (+ points) in a binary 
system with an initial orbit of P = 8 d and e = 0.70. The letters 
represent various times in the evolution of the binary system (see 
text for details). Each point represents an iteration of the BSE 
algorithm except in regions of high density, i.e. rapid evolution. 



0.2 



0.4 



Figure 2. The evolution of the spins of the primary (* points) 
and secondary (+ points) stars, and the orbital spin (x points), as 
a function of the eccentricity of the orbit for the example of Algol 
evolution, shown until the end of the first phase of mass transfer. 
The points at zero eccentricity have been slightly displaced so 
that the co-rotation of the stars with the orbit can be clearly 
seen. 



When common-envelope evolution leaves a surviving bi- 
nary system we continue the evolution with the new orbital 
parameters. If at any point the stars merge we follow the 
evolution of the resulting single star. If either star explodes 
Eis a supernova and the resultant velocity kick disrupts the 
binary we evolve the system forward as two single stars that 
don't interact. 

The CPU time required by the BSE algorithm to evolve 
500 000 binaries up to the age of the Galaxy is approximately 
20 hr, or 0.144 s per binary, on a Sun SparcUltralO worksta- 
tion (containing a 300 MHz processor). 



3 SOME EXAMPLES AND COMPARISONS 

To illustrate the BSE algorithm we present some example 
evolutionary scenarios, concentrating in particular on Al- 
gol and cataclysmic variable evolution. These examples also 
serve to demonstrate the sensitivity of binary evolution to 
the choice of model parameters where, unless otherwise spec- 
ified, the default values given in Table ^ are used. In addi- 
tion, we take a selection of evolution scenarios previously 
outlined by other authors and we compare the results of 
their models with those of the BSE algorithm. 

3.1 Algol Evolution 

Consider a system in which the initial masses are 2.9 M0 
and 0.9 M0 in an 8 d orbit with e = 0.7. In this example 
Z = 0.02 and ace = 3. Figure |^ shows the evolution of 
the binary stars in the Hertzsprung-Russell diagram. On 
the MS tidal friction arising from convective damping of 



the tide raised on the 0.9 Mq secondary acts to circular- 
ize the system so that at a time T = 413 Myr, when the 
2.9 Mq primary reaches the end of its MS lifetime, the ec- 
centricity has fallen to 0.28 and P to 3 d. The radius of 
the primary then increases rapidly as it evolves across the 
HG until Ri = 6.2 when it fills its Roche-lobe. By then 
P = 2.7 d and the orbit has circularized (position a in Fig- 
ure ^). Mass transfer proceeds on a thermal timescale, with 
the secondary accreting 80% of the transferred mass (lim- 
ited by its thermal timescale), until the primary reaches the 
end of the HG with Mi = 0.53 Mq, Mi = 2.74 M© and 
P = 17 d. During this phase the equilibrium radius of the 
primary has exceeded its Roche-lobe by as much as a factor 
of 2. Mass transfer continues as the primary begins to as- 
cend the GB, but now on a nuclear timescale which allows 
the secondary to accrete all the transferred mass, until the 
primary's envelope mass becomes so small that it shrinks in- 
side its Roche-lobe (position b) . Figure ^ shows the intrinsic 
spin of the stars and the spin of the orbit as a function of 
the eccentricity, up to this point. We see clearly that co- 
rotation is achieved as the orbit circularizes. When RLOF 
ends Ml = 0.42 M© , M2 = 2.85 Mq and P = 30 d. While the 
primary was transferring mass on the HG the mass-ratio of 
the system inverted, gi < 1, so that the more evolved com- 
ponent is now the least massive. This is the well known Algol 
Paradox (Hoyle 1955; Crawford 1955). Many semi-detached 
Algol systems have since been observed and their parame- 
ters determined (e.g. Popper 1980). These include RY Gem: 
Ml = 2.6 Mq, M2 = 0.6 Mq and P = 9.3d, and TT Hya: 
Afi = 2.6 Mq, M2 = 0.7 Mq and P = 6.95 d; both of which 
could easily be explained by this example. Figure shows 
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t3 




M/Mg 

Figure 3. The evolution of stellar mass witli binary period for 
the example of Algol evolution, until the end of the first phase 
of mass transfer. The primary mass decreases from 2.9 Mq to 
0.42 M0 (* points) and the secondary mass increases from 0.9 Mq 
to 2.85 Mq (+ points) while the period initially decreases owing 
to circularization of the orbit and then increases as a result of the 
nearly conservative mass transfer. The observed parameters for 
the Algol systems RY Gem (solid triangles) and TT Hya (solid 
squares) are also shown. 



the evolution of the mass of the stars with the orbital period 
up to the end of the Algol phase of evolution. 

Shortly after the mass transfer phase has ended the GB 
primary loses all of its envelope leaving a 0.42 Mq HeMS 
star with a 2.85 Mq MS companion which would be a blue 
straggler in a star cluster (position c). The naked helium 
star evolves to the HeGB and then to a 0.417 Mq COWD 
after losing its helium envelope (position d). The MS star, 
now the primary, evolves to the GB and fills its Roche-lobe 
(position e) when the period has reduced to 18 d as a result 
of the transfer of angular momentum from the orbit to the 
spin of the star. As qi — 6.83 > gcrit the mass transfer 
is dynamical, a CE forms and leaves a binary containing 
the COWD and the 0.416 Mq He core of the primary in an 
orbit of 0.04 d. The HeMS star evolves to fiU its Roche-lobe 
(position f ) which has shrunk by gravitational radiation, and 
transfers 0.1 Mq of He-rich material to the COWD before 
the system reaches a contact state at P = 1 min (position g). 
The two stars merge to a single HeGB star. 

3.2 A Cataclysmic Variable 

Next we use a binary with initially Mi = 6.0 Mq and 
M2 = 1.3Mq, P = 630 d, and a circular orbit, e = 0.0, 
to give a quick illustration of CV formation and evolution. 
In this example Z — 0.02, cice ~ 1 and e = 0.001. The rela- 
tively large orbital period gives the 6.0 Mq primary enough 
room to evolve on to the EAGB before filling its Roche- 
lobe 78 Myr after the evolution began. At this point the 
tides raised on the convective envelope of the primary have 



xi 



10 * 



0.5 



M/M 







Figure 4. The evolution of stellar radius (-1- points for the pri- 
mary and * points for the COWD secondary) with mass during 
a typical phase of CV evolution. The numbers indicate the evo- 
lution state of each star. Dotted lines are the Roche-lobe radii 
for each star. Note that the COWD accretes little material as its 
Roche-lobe radius shrinks and then expands. 



moderately decreased the orbital period, offset slightly by 
orbital changes owing to 0.1 Mq lost from the primary, none 
of which has been accreted by the MS companion. Because 
gi — 4.5 > Qcrit mass transfer is dynamical and a common- 
envelope forms. The He core of the primary and the MS sec- 
ondary spiral together until their orbital period is reduced 
to 0.63 d at which point the envelope has been driven off. 

The primary is now a HeHG star of mass 1.5 Mq and 
very soon evolves to fill its Roche-lobe again, resulting in 
another CE system. Once again the cores do not coalesce 
before the envelope is removed and a binary consisting of a 
0.94 Mq COWD and a 1.3 Mq MS star in a circular orbit 
of 0.45 d emerges. Gravitational radiation acts to reduce the 
orbital period until at T = 303 Myr the MS primary fills its 
Roche-lobe and a CV is formed. Mass transfer is steady and 
proceeds on a nuclear timescale, driven in part by gravita- 
tional radiation and, when the primary mass is reduced be- 
low 1.25 Mq and it develops a convective envelope, magnetic 
braking and tidal friction. As the primary evolves across the 
HG it transfers mass on a thermal timescale, now driven by 
the evolution. Throughout RLOF the mass transfer rate re- 
mains below 10~^ Mq yr"'^, low enough that nova explosions 
at the surface of the white dwarf have allowed only 0.1% of 
the transferred material to be accreted. RLOF ends with the 
primary on the GB just before it loses the remainder of its 
envelope becoming a low-mass He WD of Afi — 0.158 Mq. 
The orbital period is now 4.5 d and a gravitational radiation 
timescale of lO^** yr means that the double-degenerate sys- 
tem cannot coalesce within the age of the Galaxy. Figure ^ 
shows the radius of each star and the corresponding Roche- 
lobe radius, as a function of the stellar mass, during the CV 
phase of evolution. 
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3.3 Sensitivity to Model Parameters 

To give an idea of how sensitive tiie final state of the system 
is to changes in the initial binary parameters and the physi- 
cal parameters that govern the evolution, we first reconsider 
the Algol example described above. If the WD-HeMS sys- 
tem emerged from the CE with a separation increased by 
only a few percent then the COWD would accrete 0.15 M© 
of He-rich material and explode as an ELD SNIa before the 
system evolved into contact. Such an increase in separation 
could easily be achieved by a less than 10% change to qce, 
well within its uncertainty, or small changes in the initial 
period and/or eccentricity of the system. 

There is no doubt that major changes to the outcome 
would result if we neglected to model tidal evolution for this 
example but to do so for a system that is initially eccen- 
tric with a short period would not be physically correct. 
However, for the purpose of comparison we can examine the 
evolution of a system that begins with the same component 
masses and semi-latus rectum but in a circular orbit, i.e. 
P = 2.9 d. Apart from the initial circularization phase, this 
example evolved with the tidal evolution activated follows 
a virtually identical path to that of t he abo ve example (a 
result we will discuss further in Section 4.2.3). If tidal evolu- 
tion is not activated, and therefore any spin-orbit coupling 
is neglected, the orbit experiences a greater degree of ex- 
pansion during RLOF so that P = 25 d at the end of the 
Algol phase. The wider orbit translates to a lower mass- 
transfer rate during this phase which means that more mass 
is accreted by the companion. The subsequent evolution still 
leads to a CE when the original secondary has reached the 
GB but in this case a lower value of Qce is required if a 
COWD-HeMS binary is to emerge with P = 0.04 d. 

If the CV example presented above had a common- 
envelope efficiency of Oce = 3 less energy is required to drive 
off the CE and the WD-MS binary forms with a wider or- 
bital period of 3.0 d. This allows the MS primary to evolve 
to the GB before it fills its Roche-lobe so that dynamical 
mass transfer leads to another CE phase in which the WD 
and the giant core coalesce to form an EAGB star. A CV is 
not formed at all. If the metallicity of the component stars 
is taken to be Z = 0.001, rather than Z = 0.02, the system 
coalesces after the first instance of CE evolution and forms 
an EAGB star. This is also the case if the tidal evolution is 
not followed. 



3.4 Comparison with Other Authors 

The Algol example described by Tout et al. (1997) begins 
with initial masses 3.2 and 2.0 M© and a circular 5d orbit. 
The primary fills its Roche-lobe on the HG and mass transfer 
proceeds until it has lost its envelope on the GB, at which 
point it is a 0.318 M© naked helium star in a 170 d orbit 
with a 4.65 Mq MS companion. The MS star evolves to the 
AGB, by which time the secondary has become a COWD, 
before filling its Roche-lobe which results in CE evolution 
and the production of a double-degenerate WD binary in 
a 35min orbit. Gravitational radiation decreases the size of 
the orbit so that steady mass transfer from the less massive 
WD begins and an AM CVn system is formed. The same 
binary evolved with the BSE algorithm also forms an Algol 
system but the first phase of RLOF ends when the primary 



is in the CHeB stage and loses its envelope to become a 
0.49 Mq HeMS star with a 4.7 Mq MS companion. Tidal 
friction within the system has restrained the growth of the 
separation during mass transfer so that the orbital period 
at this stage is only 110 d. As a result the new primary fills 
its Roche-lobe on the GB and the CE evolution that follows 
leads to coalescence and the formation of a single CHeB star. 
An AM CVn system is not produced from this set of initial 
conditions when tides are included. 

Tout et al. (1997) describe a scenario for the formation 
of a CV with initial masses of 7.0 Mq and 1.0 Mq in a cir- 
cular orbit of P = 1200 d. In this example Z = 0.02 and 
Qce = 1. They find that the 7.0 Mq primary fills its Roche- 
lobe on the AGB when its mass has fallen to 5.36 Mq. The 
ensuing CE produces a 1.014 Mq COWD (they did not dis- 
tinguish ONeWDs) in a 35 hr orbit with the 1.0 Mq MS star. 
The MS star then fills its Roche-lobe, as magnetic braking 
reduces orbital separation, resulting in steady mass transfer 
with nova explosions on the WD surface. Evolution with the 
BSE algorithm leads to a similar outcome but with differ- 
ences in the states of the stars and the orbit. The 7.0 Mq 
primary fills its Roche-lobe earlier on the AGB when its mass 
has only fallen to 6.8 Mq and the CE produces a 1.183 Mq 
ONeWD in a 18 hr orbit with the 1.0 Mq MS star. Once 
again a CV with nova outbursts is formed when the MS pri- 
mary fills its Roche-lobe. In this case the differences between 
Tout et al. (1997) and the BSE algorithm are primarily due 
to improvements in the single star evolution model. 

As a final example we consider the formation of the 
double neutron star binary described by Portegies-Zwart 
& Verbunt (1996). In their model stellar evolution is per- 
formed using a simple prescription based on models with 
Z — 0.02. They assume no kick in a supernova explosion, 
and their treatment of CE evolution is similar to that of Iben 



Livio (1993) discussed in Section 2.7.1. Tidal evolution is 



not included in their model but eccentric orbits are instan- 
taneously circularized when the radius of either star exceeds 
1/5 of the separation at periastron. The initial binary has 
masses 13.1 and 9.8 Mq separated by 138 Rq in a circular 
orbit. Mass transfer, which begins when the primary has 
evolved to the HG, becomes dynamical when the primary 
reaches the GB leading, via CE evolution, to the formation 
of a binary composed of a 3.7 Mq He star and an 18.7 Mq 
MS companion separated by 411 Rq. The He star then un- 
dergoes a supernova which leaves a 1.34 Mq NS remnant 
and widens the orbital separation to 463 Rq . Another phase 
of dynamical mass transfer, when the 18.7 Mq star fills its 
Roche-lobe, leads to CE from which the NS emerges with a 
4.6 Mq He star companion, separated by 1.6 Rq. When the 
new He star explodes as a supernova a double NS binary 
forms with P — 0.5 d. 

This example evolved with the BSE algorithm forms a 
binary consisting of a 3.1 Mq HeMS star separated from an 
11.7 Mq MS companion by 60 Rq after the first CE phase. 
The separation in this instance is much smaller because not 
as much mass is transferred and because tides cause the 
spin of each star to synchronize with the binary period so 
that the orbit does not expand during mass transfer. Next 
the HeMS star evolves to the HeGB before a supernova ex- 
plosion leaves a 1.33 Mq NS and widens the separation to 
67 Rq. The second phase of CE evolution then leads to co- 
alescence of the NS and the 2.63 Mq He core which cre- 
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ates an unstable Thorne-Zytkow object. However, if the ini- 
tial binary consists of masses 13.6 and 7.7 separated by 
100 R© , the evolution proceeds along a similar path to that 
described by Portegies-Zwart & Verbunt (1996), forming a 
double NS binary with P = 1 d. This system coalesces in 
300 Myr through gravitational radiation. If a velocity kick 
at supernova is included in the model then the final state of 
the system depends on the magnitude of the kick and the 
relative positions of the two stars when the explosion occurs. 



4 THE EFFECT OF TIDES ON BINARY 
EVOLUTION 

The examples in the previous section have shown that the 
evolution of a binary is sensitive to the inclusion of tidal evo- 
lution in the binary algorithm. This is true even for initially 
circular orbits. It is instructive to investigate and quantify 
the systematic effect of tides within a population of bina- 
ries. We do this using the rapid evolution code to generate 
a series of large binary populations from which we calculate 
the formation rate of interesting individual species, such as 
X-ray binaries, double-degenerates and symbiotic stars, and 
events such as type la SNe. The results for populations with 
and without tidal friction can be compared. In addition it 
is useful to investigate the influence of various uncertain 
parameters in the binary algorithm, such as the common- 
envelope efficiency parameter. Comparison with observa- 
tions and work of others hopefully constrains the value of 
some of these parameters. The statistics of theoretical binary 
populations have been presented and discussed by many au- 
thors in the past. Examples include Han (1998): the forma- 
tion of double-degenerates, type la SNe and CVs; Portegies 
Zwart & Yungelson (1998): the formation of binary neu- 
tron stars; de Kool (1992) and Politano (1996): formation 
of CVs; Yungelson et al. (1995): the formation of symbiotic 
stars with a WD secondary; Iben, Tutukov & Yungelson 
(1995): low- mass X-ray binaries in the Galactic disk; Han et 
al. (1995): the formation of CVs, Algols, double-degenerates 
and symbiotic stars; and Pols & Marinus (1994); the produc- 
tion of blue stragglers. These studies are all flawed by the 
fact that no attempt was made to model the tidal circular- 
ization and synchronization of the binary orbit as it evolves, 
and in particular, the interaction of the intrinsic spin of the 
stars with the orbital period. 

4.1 Method 

Our aim is to evolve a population of binaries according to 
chosen distributions of primary mass, secondary mass and 
orbital separation, in conjunction with a realistic birth rate 
function and, from this population, to calculate birth rates 
and expected numbers in the Galaxy for various individual 
binary populations and interesting events. 

We flrst set-up a grid of initial binary parameters Mi, 
M2 and a within the limits: 



S\nX = 



Tlx 



[In Xmax — In A'min 



(93) 



Ml 


0.8 - 


80 Mq 


M2 


0.1 - 


80 M© 


a 


3.0 - 


^ 10* R© 



with the nx grid points of parameter X logarithmically 
spaced. 



For each set of initial parameters we evolve the binary sys- 
tem to an age of 15 Gyr, or until it is destroyed. Each phase 
of the evolution is followed in detail according to the BSE al- 
gorithm. The lower limit on the primary mass is determined 
by the lowest mass star that will change appreciably in size 
within 15 Gyr and systems with M2 > Mi are not evolved. 
If a particular binary system j evolves through a phase that 
is to be identifled with a certain individual binary popula- 
tion i, such as cataclysmic variables, then the system makes 
a contribution 

Srj = 5$ (In Mij) ifi (In Msj) * (In aj) 5 In A/i5 In MiS In a(94) 

to the rate r; at which that particular population is born. 
This rate depends on the star formation rate S, the primary 
mass distribution "I>(lnMi), the secondary mass distribu- 
tion (p (lnM2), and the separation distribution ^ (In a). The 
number of population i in the Galaxy at any time T is then 
given by the sum of 5rj x Aiji for all systems j that lived 
for a time Atji as a member of that population, where we 
assume that the interval began before T and that T is the 
endpoint of the binary evolution. 



4-1.1 The Mass Distributions 

Our primary mass distribution is the IMF of Kroupa, Tout & 
Gilmore (1993), derived from the stellar distribution towards 
both Galactic poles as well as the distribution of stars within 
5.2 pc of the Sun, 







m < mo 



, N _ I aim mo < m < 0.5 
^^"^'^S a^m-^-^ 0.5<m<1.0 
a2m~^''^ 1.0 < m < 00 



(95) 



where ^ (m) dm is the probability that a star has a mass, 
expressed in solar units, between m and m -I- dm. The dis- 
tribution is normalized according to 



^(m) dm - 



(96) 



so that, for mo = 0.1, ai = 0.29056 and 02 = 0.15571. Then 

<l>(lnMi) = Mie(Mi) . (97) 

The percentage of stars with mass greater than 80 M© is less 
than 0.005 and hence our upper mass limit for the grid. 

If the component masses are to be chosen independently 
from the IMF then the secondary mass distribution is 



V9(lnM2) =M2C(M2) 



(98) 



However there is observational evidence (Eggleton, Fitchett 
& Tout 1989; Mazeh et al. 1992; Goldberg & Mazeh 1994) 
to support correlated masses, i.e. 



(^(lnM2) = = 92 , 



(99) 



which corresponds to a uniform distribution of the mass- 
ratio 52, for < 52 < 1. 

4.1.2 The Separation Distribution 
We take this to be 



© 2000 RAS, MNRAS 000, 



20 J. R. Hurley, C. A. Tout and 0. R. Pols 



*(lna) 



constant, 



(100) 



between the limits 3 and 10'' R©. Normalization gives k — 
0.12328. This distribution is contrary to the findings of 
Eggleton, Fitchett & Tout (1989) but it facifitates compar- 
ison with other models (e.g. Yungelson, Livio & Tutukov 
1997; Han 1998) without upsetting the results. We chose 
the upper limit with a view to including all separations for 
which a binary system is likely to experience some form of 
mass exchange interaction within the lifetime of the Galaxy 
(Yungelson et al. 1995) . This was not achieved for all classes 
of binaries (see Section 4.2.2). The choice of lower limit is 
rather simplistic and a dependence on binary mass might be 
more realistic. 



4- 1-3 Star Formation Rate 

We assume that one binary with Afi > 0.8 Mq is born in 
the Galaxy per year. Therefore 

/• All — ^ 

S $(m) dm = 1, (101) 

J J\/l=0.8 

which gives S = 7.6085 yr~^. We fix this rate over the life- 
time of the Galaxy. It is in rough agreement with the birth 
rate of WDs in the Galaxy, Xwd ~ 2 x 10~'^ pc~^yr~^ 
(Phillips 1989), when we note that only stars with A/ > 
0.8 M0 can evolve to white dwarfs in the age of the Galaxy, 
that all stars are in binaries and assume an efi'ective Galactic 
volume of Vgai = 5 x lO''^ pc^. The same assumptions regard- 
ing the star formation rate have been made previously (e.g. 
Iben & Tutukov 1984; Han 1998) facilitating comparison 
with the results of these studies. 

To estimate an effective Galactic volume we use val- 
ues for the Galactic disk of I/v.disk = 1.2 X 10'° L© and 
(M/L)^j^j^ = 5 for the visual luminosity and mass-to-light 
ratio (Binney & Tremaine 1987) to give AfdiA = 6x10'° Mo 
for the mass in visible stars. We combine this with the lo- 



cal mass density of stars 
Gilmore 1989), to 



P© 



ive Fgai = 6 X 10" pc^ 



0.1 Mq pc"'^ (Kuijken & 



4-1-4 Binary Population Models 

We construct a variety of models, each with slightly differ- 
ent assumptions for the initial conditions of the population 
or the parameters that govern the evolution. Our standard 
is Model A in which the secondary mass is chosen according 
to a uniform distribution of the mass-ratio (see eq. p9[)and 
the metallicity of the stars is Z — 0.02. All binaries are 
initially circular, with each star having ZAMS rotation as 
given by eq. (108) of Papl, and tidal synchronization is fol- 
lowed in detail. The common-envelope efficiency parame- 
ter is set at acE = 3 so that the treatment of common- 
envelope evolution in this model is similar to that of Iben & 
Livio (1993, although the dependence on core and envelope 
mass in the two treatments remains somewhat different). 
Model B differs from Model A by not including tidal evolu- 
tion. For Model C the tides are taken into account but we 
use QfcE = 1. Model D is the same as Model A except that 
the secondary mass is chosen independently from the same 
IMF as the primary. 

The assumption that all stars within the population are 



Table 4. Model parameters. 



TIDES acE 



A/2 



A 
B 
C 
D 
E 
F 
G 



ON 
OFF 
ON 
ON 
ON 
ON 
PZV 



3.0 
3.0 
1.0 
3.0 
3.0 
3.0 
3.0 



0.0 
0.0 
0.0 

0.0 
0.0 

fie) 
/(e) 



eq. 
eq. 
eq. 
eq. 
eq. 
eq. 
eq. 



(||) 
(i) 
(i) 



0.0200 
0.0200 
0.0200 
0.0200 
0.0001 
0.0200 
0.0200 



born with the same composition is somewhat naive: nucle- 
osynthesis in successive generations of stars enriches the gas 
from which they form as the Galaxy evolves. Currently it is 
uncertain whether a definite age-metallicity relation exists 
for stars in the Galactic disk but it is clear is that there is 
a large scatter in metallicity at any given age (McWilliam 
1997). From the data of Edvardsson et al. (1993) the age- 
metallicity relation for the solar neighbourhood can be rep- 
resented in terms of the relative iron abundance by 



[Fe/H] = 0.35a; - 0.4r - (1.73 + 0.352:) t 



(102) 



where r = t/T^ai is the fractional time since formation of 
the star in terms of the current Galactic age. A star born 
now has r = 0. The scatter in the data is reproduced by a 
random variable x uniform in the range — 1 < a; < 1. The 
metallicity can be obtained from 



[Fe/H] ~ln(^/Ze 



(103) 



with Zq — 0.02. This relation shows that solar composition 
is within the limits of the scatter for the most recent 80% of 
the Galaxy's life but that for older populations the metal- 
licity drops sharply to much lower values. While it would 
be interesting to construct a model with a Z{t) dependence 
this is beyond the scope of this work. To represent an older 
population we simply take the metallicity of the stars in 
Model E to be ^ = 0.0001. 

Finally to investigate the changes produced when tidal 
circularization, as well as synchronization, is acting on the 
population, our Model F differs from Model A by allowing 
binaries to form in eccentric orbits. We choose the initial 
eccentricity from a thermal distribution (Heggie 1975) 



/(e) = 2e 



(104) 



between the limits and 1. Model G also allows eccentric 
orbits but the process of circularization is according to the 
model of Portegies Zwart & Verbunt (1996, PZV), as de- 
scribed in Section ^. 

The characteristics of each model are summarised in 
Table ^. All models assume the default values given in Ta- 
ble ^ unless specified otherwise. This means that a velocity 
kick is imp arted to a NS or BH at birth (as described in 
Section |2.5| ) and the Eddington limit is not imposed during 
mass transfer in any of these models (see Section 2.6.6). We 
set the wind accretion parameter to /3w = 0.5 in all cases 
and the stars are not assumed to be in co-rotation with the 
orbit on the ZAMS. 

For each model we use 100 grid points in each dimension 
so that, recalling that only systems with ^2 < 1 need be 
evolved, a total of 6.54 x 10^ different binaries are evolved 
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in each case. We find that a finer grid spacing does not alter 
the results within the tolerance set by the quoted errors. For 
each system we slightly displace the initial parameters from 
their grid positions so that a smooth distribution in mass 
and separation is achieved. This is done uniformly about 
each grid point, after selection of a random variable, with 
no overlap between parameters attached to adjoining grid 
points. 



4.2 Results 

Table ^ lists the formation rates, as number per year in 
the Galactic disk, of various interesting stars, binaries and 
events. We calculate these assuming an age for the Galac- 
tic disk of 12 Gyr (Eggleton, Fitchett & Tout 1989; Phelps 
1997; Knox, Hawkins & Hambly 1999) and that all stars are 
born in binaries. Also listed is an estimate of the percentage 
error in the formation rate of the individual populations cal- 
culated from the standard deviation produced by repeated 
runs of Model A, each with a different seed for the ran- 
dom number generator. Errors associated with systems that 
involve the formation of a NS or BH are generally larger 
owing to the velocity kick distribution. The eccentricity dis- 
tribution used in Models F and G introduces an additional 
degree of uncertainty in the results, however repeated runs 
of Model F show that the percentage error only exceeds that 
quoted in Table ^ for a few cases, and never by more than 
a factor of two. 



4-2.1 Definitions of Binary Class 

In our models a blue straggler star (BSS) is a MS star that, 
by accreting mass or by merging, appears older than the 
standard MS lifetime applicable to its mass. This phase of 
evolution lasts until the star moves off the MS. Cataclysmic 
variables (CVs) have a WD secondary with a non-degenerate 
Roche-lobe filling companion. These we divide into sub- 
categories, classical CV (CV class) if fei < 1, GK Persei 
systems (GK Per, e.g. Crampton, Cowley & Fisher 1986) if 
fci — 2, symbiotic-like binaries (CV Symb) if 3 < fci < 6 and 
subdwarf B binaries (sdB) if 7 < fci < 9. Binaries in which 
a MS secondary is accreting from a Roche-lobe filling com- 
panion are termed Algol. If the mass-ratio of the primary, 
qi = M1/M2, is greater than unity then the system is pre- 
Algol, if the primary is a MS star and gi < 1 then this is a 
MS Algol, otherwise the system is a cold Algol, if the sec- 
ondary has M2 < 1.25 M0, or a hot Algol, if M2 > 1.25 M©. 

X-ray binaries have a NS or a BH secondary accreting 
material from either the stellar wind of the companion or by 
RLOF, with an accretion luminosity Lx > L©. The primary 
star may be of any type. If the material is accreted via an 
accretion disk the accretion luminosity is given by 



log 



1.62 + 1.07 (P/hr) 
2.22 + 1.07 (P/hr) 



NS secondary 
BH secondary 



,(106) 



GM2M2 
2R2 



(105) 



If the primary is a MS star with Mi < 2 Mq the system is 
a low-mass X-ray binary (LMXB), if the primary is a WD 
then it is a white dwarf X-ray binary (WDXB), otherwise 
it is a massive X-ray binary (MXRB). These sources are 
divided into transient (t) or persistent (p) according to the 
criteria of van Paradijs (1996): for a critical luminosity 



if Lx > ix.crit the X-ray source is persistent, otherwise it is 
a soft X -ray t ransient (SXT) analogous to a dwarf nova (see 
Section 2.6.6). SXT outbursts are caused by an instability 



in the accretion disk, the temperature of which is affected 
by X-ray heating. The main effect of this X-r ay heating, 
taken into account in the calculation of eq. (106), is to make 



the flow of material through the accretion disk stable to 
substantially lower mass transfer rates than for dwarf novae 
systems. 

Symbiotic stars are classed as systems with secondaries 
of k2 < 12 that accrete material from the stellar wind of a 
giant primary at a rate high enough to produce an accre- 
tion luminosity that exceeds either 10 L© (Yungelson et al. 
1995) or 1% of th e primary luminosity - whichever is the 
smaller. Eq. ( [l05| ) is used to calculate the accretion lumi- 
nosity. D-type symbiotics (D-Symb) are long-period systems 
containing a Mira-like cool star primary (fci = 6). They are 
characterized by extensive circumbinary dust shells. Shorter 
period symbiotics with a normal giant primary (fci < 6) are 
S-type (S-Symb). MS stars accreting from a Roche- lobe fill- 
ing giant could also appear to be symbiotic stars (Kenyon 
1986) but here we include them in the cold Algol popula- 
tion. Miscellaneous RLOF systems with a non-degenerate 
primary (fci < 10) are nHe MSG if the secondary is a he- 
lium star (7 < fc2 < 9) and gnt MSG if the secondary is 
giant-like (2 < fc2 < 6). 

Double-degenerate (DD) binaries that consist of two 
WDs are WDWD DD while those containing a WD and 
either a NS or BH are WDNS DD. If the system is com- 
posed of a NS-NS, NS-BH or BH-BH pair it is a NSNS DD. 
DD systems with a Roche-lobe filhng WD (AM CVn sys- 
tems) are He DDRch, if the the primary is a He WD, or CO 
DDRch otherwise. NS-NS binaries that evolve into contact 
and coalesce are NS DDRch. As a result of binary interac- 
tion it is possible to form WDs of mass less than the lowest 
mass WD that can be formed from single star evolution in 
the lifetime of the Galaxy, M ^ 0.5 M0. These low-mass 
white dwarfs (LMWDs) are recorded as either He LMWD 
or CO LMWD. 

We also count various types of supernovae. HeWDs that 
explode when their mass exceeds 0.7 Mq are recorded as He 
SNIa, as are He WD mergers. COWDs that explode as a pos- 
sible type la SN in an Edge-Lit Detonation, once 0.15 M© 
of He-rich material has been accreted by the WD, are ELD 
SNIa, while GOWDs that explode because their mass ex- 
ceeds the Ghandrasekhar limit are CO SNIa. Supernovae 
that leave no remnant from a primary with fci < 9 are SNIIa. 
AIG represents the accretion induced collapse of a Ghan- 
drasekhar mass ONeWD to a NS. Supernovae that produce 
a NS or BH from a primary with 8 < fci < 9 are SNIb/c. 
All others are normal type II SNe, or SNII. 

Supernovae are effectively instantaneous events so there 
is no record of how long they last, whereas the remainder 
of the populations are systems and must exist for a non- 
zero amount of time to contribute to the birth rate. We also 
count double NS binaries that enter RLOF (NS DDRch) as 
events because they always coalesce immediately as possible 
7-ray bursts. Each binary may contribute to the rate of more 
than one population if it evolves through a series of phases 
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Table 5. Formation rates per year in the Galactic disk of various events and systems. 



Model 




A 










B 








C 








D 








E 








F 








G 






%err 




BSS 


1.295 


X 


10" 


-1 


1 


.138 


X 


10" 


-1 


1.248 


X 


10" 


-1 


5.049 


X 


10" 


-2 


1.839 


X 


10" 


-1 


1.181 


X 


10" 


-1 


1.367 


X 


10" 


-1 


0.31 


CV class 


1.987 


X 


10" 


-2 


1 


.292 


X 


10" 


- 2 


2.120 


X 


10" 


- 2 


1.820 


X 


10" 


-2 


3.080 


X 


10" 


-2 


1.889 


X 


10" 


- 2 


3.973 


X 


10" 


- 2 


3.65 


GK Per 


1.398 


X 


10" 


-2 


1 


.165 


X 


10" 


- 2 


1.671 


X 


10" 


- 2 


1.679 


X 


10" 


-3 


3.709 


X 


10" 


-2 


1.374 


X 


10" 


- 2 


3.460 


X 


10" 


- 3 


0.93 


CV Symb 


7.858 


X 


10" 


-4 


1 


.423 


X 


10" 


-4 


2.450 


X 


10" 


- 3 


1.898 


X 


10" 


-4 


1.027 


X 


10" 


-2 


8.527 


X 


10" 


-4 


4.440 


X 


10" 


-4 


2.73 


sdB 


1.135 


X 


10" 


-2 


1 


.398 


X 


10" 


-2 


2.189 


X 


10" 


-3 


1.090 


X 


10" 


-3 


1.729 


X 


10" 


-2 


1.047 


X 


10" 


~2 


1.125 


X 


10" 


"-^ 


0.49 


pre Algol 


1.248 


X 


10~ 


_ 1 


1 


.040 


X 


10" 


_ 1 


1.244 


X 


10" 


_ 1 


8.807 


X 


10" 


_2 


1.861 


X 


10" 


_ 1 


1.112 


X 


10" 


_ 1 


2.082 


X 


10" 


_ 1 


0.07 


MS Algol 


2.861 


X 


10' 


_2 


2 


.981 


X 


10" 


_ 2 


2.864 


X 


10" 


_ 2 


1.558 


X 


10" 


_2 


2.784 


X 


10" 


_2 


2.485 


X 


10" 


_ 2 


1.077 


X 


10" 


_ 2 


0.32 


cold Algol 


1.313 


X 


10" 


-2 


1 


.145 


X 


10" 


_ 2 


1.293 


X 


10" 


_ 2 


7.007 


X 


10" 


-3 


3.718 


X 


10" 


-2 


1.222 


X 


10" 


_ 2 


2.255 


X 


10" 


_ 3 


1.11 


hot Algol 


4.910 


X 


10" 


-2 


4, 


.593 


X 


10" 


-2 


4.458 


X 


10" 


-2 


1.018 


X 


10" 


-2 


9.049 


X 


10" 


-2 


4.430 


X 


10" 


-2 


2.780 


X 


10" 


-2 


0.18 


NS LMXBp 


3.487 


X 


10" 


-6 


1 


.362 


X 


10" 


_ 7 


2.219 


X 


10" 


_ 7 


5.481 


X 


10" 


-6 


1.010 


X 


10" 




1.071 


X 


10" 


-Q 


4.906 


X 


10" 


-Q 


37.56 


BH LMXBp 


3.601 


X 


10" 


-6 


1 


.696 


X 


10" 


-6 


1.561 


X 


10" 


-6 


3.686 


X 


10" 


~6 


2.563 


X 


10" 


-5 


3.865 


X 


10" 


-6 


5.349 


X 


10" 




25.36 


NS MXRBp 


7.481 


X 


10" 


-4 


5 


.014 


X 


10" 


-4 


6.307 


X 


10" 


-4 


1.434 


X 


10" 


~5 


6.096 


X 


10" 


-4 


7.214 


X 


10" 


-4 


5.164 


X 


10" 


-4 


1.37 


BH MXRBp 


1.057 


X 


10" 


-4 


1 


.612 


X 


10" 


-5 


1.092 


X 


10" 


-4 


2.043 


X 


10" 


-6 


2.425 


X 


10" 


-4 


1.085 


X 


10" 


-4 


1.132 


X 


10" 


-4 


3.77 


IN O VV 


1 639 




10" 


-3 


2 


028 




10" 


-3 


1 686 




10" 


~4 


8 393 




10" 


-5 


3 563 




10" 


-3 


1 426 




10" 


-3 


1 347 




10" 


-3 


55 


BH WDXBp 


2.764 


X 


10" 


-4 


3, 


.513 


X 


10" 


-4 


2.885 


X 


10" 


~5 


4.987 


X 


10" 


-6 


4.951 


X 


10" 


-4 


2.400 


X 


10" 


-4 


2.085 


X 


10" 


-4 


7.62 


NS LMXBt 


2.359 


X 


10" 


-5 


5, 


.640 


X 


10" 


-6 


1.253 


X 


10" 


-6 


2.537 


X 


10" 


-5 


4.198 


X 


10" 


-5 


9.087 


X 


10" 


-6 


3.654 


X 


10" 


-5 


21.27 


BH LMXBt 


9.172 


X 


10" 


-6 


1 


.691 


X 


10" 


"6 


4.854 


X 


10" 


'6 


6.552 


X 


10" 


-6 


2.323 


X 


10" 


-5 


6.542 


X 


10" 


-6 


2.283 


X 


10" 


-5 


20.87 


NS MXRBt 


7.345 


X 


10" 


-4 


8 


.048 


X 


10" 


'4 


8.122 


X 


10" 


-4 


1.823 


X 


10" 


_5 


1.568 


X 


10" 


„3 


7.149 


X 


10" 


-4 


6.515 


X 


10" 


-4 


1.41 


BH MXRBt 


5.447 


X 


10" 


-5 


2 


.320 


X 


10" 


-5 


8.706 


X 


10" 


'5 


2.564 


X 


10" 


~6 


2.560 


X 


10" 


-4 


5.823 


X 


10" 


-5 


4.471 


X 


10" 


-5 


2.81 


NS WDXBt 


8.963 


X 


10" 


-4 


1 


.088 


X 


10" 


-3 


7.332 


X 


10" 


-5 


7.325 


X 


10" 


-5 


2.296 


X 


10" 


-3 


7.753 


X 


10" 


~4 


6.918 


X 


10" 


-4 


0.29 


BH WDXBt 


6.531 


X 


10" 


-4 


8, 


.299 


X 


10" 


-4 


8.357 


X 


10" 


-5 


1.017 


X 


10" 


-5 


1.237 


X 


10" 


-3 


5.636 


X 


10" 


~4 


5.784 


X 


10" 


-4 


1.78 


S-Syiiib 


5.353 


X 


10" 


-3 


5 


.091 


X 


10" 


-3 


5.370 


X 


10" 


-3 


2.712 


X 


10" 


-4 


6.356 


X 


10" 


-3 


4.100 


X 


10" 


-3 


4.305 


X 


10" 


-3 


1.79 


D-Symb 


4.322 


X 


10" 


-2 


4, 


.701 


X 


10" 


-2 


4.302 


X 


10" 


-2 


5.748 


X 


10" 


-3 


3.782 


X 


10" 


-2 


3.494 


X 


10" 


-2 


3.763 


X 


10" 


-2 


0.23 


nHc MSG 


6.441 


X 


10" 


-3 


1 


.929 


X 


10" 


-3 


2.160 


X 


10" 


-3 


1.657 


X 


10" 


-3 


1.649 


X 


10" 


-2 


5.591 


X 


10" 


-3 


5.059 


X 


10" 


-3 


0.94 


gnt MSG 


3.366 


X 


10" 


-3 


3 


.639 


X 


10" 


-3 


3.432 


X 


10" 


-3 


1.093 


X 


10" 


-3 


4.809 


X 


10" 


-3 


3.075 


X 


10" 


-3 


8.166 


X 


10" 


-4 


4.75 


WDWD DD 


1.131 


X 


10" 


~1 


1 


.229 


X 


10" 


-1 


7.572 


X 


10" 


-2 


1.334 


X 


10" 


-2 


2.290 


X 


10" 


-1 


8.631 


X 


10" 


~2 


7.902 


X 


10" 


-2 


0.16 


WDNS DD 


8.577 


X 


10" 


-4 


8, 


.982 


X 


10" 


-4 


5.265 


X 


10" 


-4 


2.088 


X 


10" 


-5 


2.842 


X 


10" 


-3 


7.986 


X 


10" 


-4 


7.537 


X 


10" 


-4 


1.10 


NSNS DD 


7.384 


X 


10" 


-5 


8, 


.192 


X 


10" 


-5 


2.785 


X 


10" 


~5 


6.972 


X 


10" 


-7 


1.604 


X 


10" 


-4 


7.587 


X 


10" 


-5 


7.295 


X 


10" 


-5 


5.11 


He DDRchj^ 


1.997 




10" 


-2 


2 


.385 




10" 


-2 


6.782 




10" 


-3 


2.344 




10" 


-3 


5.361 




10" 


-2 


1.927 




10" 


-2 


1.619 




10" 


-2 


1.21 


GO DDRch 


2.942 


X 


10" 


-3 


3, 


.598 


X 


10" 


-3 


1.549 


X 


10" 


~4 


8.558 


X 


10" 


-5 


4.959 


X 


10" 


-3 


2.538 


X 


10" 


-3 


2.144 


X 


10" 


-3 


0.43 


NS DDRch 


5.724 


X 


10" 


-5 


6, 


.826 


X 


10" 


-5 


2.381 


X 


10" 


-5 


5.096 


X 


10" 


-7 


1.062 


X 


10" 


-4 


5.884 


X 


10" 


-5 


5.240 


X 


10" 


-5 


6.12 


He LMWD 


1.671 


X 


10" 


-1 


1 


.709 


X 


10" 


-1 


9.715 


X 


10" 


-2 


1.138 


X 


10" 


-1 


2.985 


X 


10" 


-1 


1.658 


X 


10" 


-1 


1.466 


X 


10" 


-1 


0.42 


GO LMWD 


7.019 


X 


10" 


-3 


7, 


.622 


X 


10" 


-3 


1.351 


X 


10" 


-3 


1.564 


X 


10" 


-3 


1.085 


X 


10" 


-2 


6.597 


X 


10" 


-3 


3.002 


X 


10" 


-3 


1.78 


He SNIa 


2.818 


X 


10" 


-3 


8, 


.100 


X 


10" 


-3 


3.217 


X 


10" 


-4 


7.214 


X 


10" 


-4 


2.702 


X 


10" 


-2 


2.677 


X 


10" 


-3 


1.242 


X 


10" 


-3 


1.68 


ELD SNIa 


1.560 


X 


10" 


-2 


1 


.994 


X 


10" 


-2 


5.193 


X 


10" 


-3 


1.665 


X 


10" 


-3 


3.227 


X 


10" 


-2 


1.490 


X 


10" 


-2 


1.102 


X 


10" 


-2 


0.61 


GO SNIa 


2.567 


X 


10" 


-3 


2 


.627 


X 


10" 


-3 


2.109 


X 


10" 


-4 


1.125 


X 


10" 


-4 


2.720 


X 


10" 


-3 


2.442 


X 


10" 


-3 


1.890 


X 


10" 


-3 


0.59 


AIC 


3.738 


X 


10" 


-3 


3, 


.975 


X 


10" 


-3 


9.771 


X 


10" 


~4 


1.364 


X 


10" 


-4 


6.286 


X 


10" 


-3 


3.442 


X 


10" 


-3 


3.029 


X 


10" 


-3 


0.35 


SNII 


1.724 


X 


10" 


-2 


1 


.592 


X 


10" 


-2 


1.940 


X 


10" 


-2 


1.142 


X 


10" 


-2 


2.848 


X 


10" 


-2 


1.810 


X 


10" 


-2 


1.921 


X 


10" 


-2 


1.75 


SNIIa 


3.551 


X 


10" 


-4 


1 


.453 


X 


10" 


-4 


6.189 


X 


10" 




7.054 


X 


10" 


-6 


6.840 


X 


10" 


-5 


2.851 


X 


10" 


-4 


1.228 


X 


10" 


-4 


1.67 


SNIb/c 


1.403 


X 


10" 


-2 


1 


.361 


X 


10" 


-2 


1.304 


X 


10" 


-2 


4.043 


X 


10" 


-3 


2.026 


X 


10" 


-2 


1.313 


X 


10" 


~2 


1.510 


X 


10" 


- 2 


0.19 



" DDRch systems are a subset of the corresponding DDs. 



identified with the individual populations listed in Table |5| 
The expected numbers presently in the Galactic disk for 
all the individual populations that exist as systems for a 
finite time are listed in Table ^. These can be converted to 
space densities by dividing by the effective Galactic volume 
Vgai = 6 X 10^^ pc^ which has been normalized to the local 
mass density of stars. 

4-2.2 Comparison of Model Results 

Comparison of Models A and B reveals that most popu- 
lations show differences in their formation rates that are 
greater than the errors but surprisingly few exhibit signifi- 
cant differences. Figure ^ shows the parameter space for the 
formation of all CV types when the initial mass of one com- 
ponent star is fixed at 3.2 Mq. The systems that form in 
Model A only, Model B only and in both are distinguished. 



There is no special reason for the choice of 3.2 Mq as an il- 
lustrative example other than that a wide variety of systems 
are formed in this case. Systems with low-mass ZAMS com- 
panions (M2 < 0.5 M0) are CV class while GK Per systems 
are formed with high-mass ZAMS primaries (Mi > 7.5 Mq). 
The sdB systems form in two main regions of the parame- 
ter space: small orbital separation (a < 30 Rq) and large 
orbital separation (a > 500 R0). The remainder, those oc- 
cupying the middle region of the diagram, are either CV 
class, GK Per or CV Symb, with a large number of the 
systems evolving through a combination of these phases. It 
is evident from Figure a that the CV population behaves 
as expected when tidal evolution is allowed in the model. 
Wider systems are brought into an interaction distance by 
the action of tidal synchronization while closer systems, that 
formed CVs without tides, now get too close and coalesce 
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Table 6. Present number in the Galactic disk of various systems. 



Model 


A 




B 




G 




D 




E 




F 




G 




BSS 


6.999 


X 


10^ 


4.444 


X 


10'' 


6.993 


X 


10' 


4.439 


X 


10' 


8.823 


X 


10' 


8.023 


X 


10' 


6.128 


X 


10' 


CV class 


2.532 


X 


10^ 


3.594 


X 


10' 


3.554 


X 


10' 


9.908 


X 


10' 


3.995 


X 


10' 


2.550 


X 


10' 


2.565 


X 


10' 


GK Per 


7.242 


X 


10^ 


2.775 


X 


10« 


8.155 


X 


10^ 


1.184 


X 


lO'' 


2.001 


X 


10® 


6.953 


X 


10® 


9.842 


X 


10^ 


CV Symb 


7.070 


X 


IQS 


1.255 


X 


10^ 


2.503 


X 


10*5 


1.849 


X 


10^ 


4.781 


X 


10* 


7.982 


X 


10^ 


2.486 


X 


10^ 


sdB 


3.298 


X 


10-* 


6.635 


X 


10* 


5.548 


X 


10^ 


3.335 


X 


10^ 


5.054 


X 


10* 


3.047 


X 


10* 


8.599 


X 


10^ 


pr6 Aigoi 


6.396 


X 


-LU 


1.262 


X 


-LU 


6.405 


X 


-LU 


3.122 


X 


1 n6 


2.411 


X 


-LU 


5.614 


X 


1 

-LU 


7.313 


X 


1 n5 


MCI Alfrnl 


6.081 


X 


1 


1.178 


X 


in' 

-LU 


6.077 


X 


1 


3.423 


X 


1 

J-U 


3.806 


X 


1 

-LU 


5.432 


X 


1 

-LU 


2.115 


X 


1 

J-U 


cold Algol 


6.784 


X 


1 

lU 


1.020 


X 


lU 


6.795 


X 


1 n6 


4.658 


X 


1 

±u 


9.008 


X 


lU 


6.007 


X 


lU 


4.516 


X 


lU 


noi rtigoi 


7.132 


X 


1 

iU 


8.216 


X 


lU 


7.148 


X 


1 n6 

lU 


3.680 


X 


lU 


6.790 


X 


lU 


6.320 


X 


lU 


5.034 


X 


lU 


NS LMXBp 


4.191 


X 


iqo 


1.112 


X 


lOi 


3.451 


X 


10° 


7.690 


X 


lOi 


4.409 


X 


10^ 


1.331 


X 


10° 


4.388 


X 


lOi 


BH LMXBp 


9.140 


X 


IQi 


4.106 


X 


10^ 


5.016 


X 


lOi 


3.001 


X 


10^ 


8.553 


X 


10^ 


6.998 


X 


lOi 


2.196 


X 


10^ 


NS MXRBp 


2.810 


X 


10^ 


2.608 


X 


10^ 


1.600 


X 


10^ 


8.178 


X 


10^ 


3.987 


X 


10^ 


2.531 


X 


10^ 


2.494 


X 


10^ 


BH MXRBp 


2.335 


X 


10^ 


2.992 


X 


10^ 


1.439 


X 


10^ 


3.333 


X 


10° 


1.876 


X 


10^ 


2.673 


X 


10^ 


1.288 


X 


10^ 


IN o VV UAop 


1.233 


X 


lU 


1.494 


X 


lU 


1.116 


X 




7.550 


X 


lU 


2.871 


X 


lu 


1.055 


X 


lU 


1.023 


X 


1 n5 
iU 


BH WDXBp 


1.554 


X 


10-* 


2.145 


X 


10* 


1.841 


X 


10^ 


3.066 


X 


10^ 


2.947 


X 


10* 


1.469 


X 


10* 


1.353 


X 


10* 


NS LMXBt 


2.854 


X 


10* 


1.871 


X 


10* 


4.131 


X 


10^ 


1.396 


X 


10^ 


5.545 


X 


10* 


2.238 


X 


10* 


1.760 


X 


10^ 


BH LMXBt 


7.240 


X 


10* 


1.085 


X 


10* 


4.260 


X 


10* 


5.072 


X 


10* 


1.343 


X 


10^ 


4.405 


X 


10* 


2.456 


X 


10^ 


MQ MXRRt 


1.785 


X 


10^ 


1.511 


X 


10^ 


5.432 


X 


10* 


8.452 


X 


10^ 


7.547 


X 


10* 


1.630 


X 


10^ 


2.962 


X 


10^ 


BH MXRBt 


4.829 


X 


10* 


7.947 


X 


10^ 


9.105 


X 


10* 


2.818 


X 


10^ 


7.353 


X 


10* 


6.671 


X 


10* 


1.432 


X 


10^ 


NS WDXBt 


2.802 


X 


10« 


3.406 


X 


10^ 


2.079 


X 


10^ 


2.031 


X 


10^ 


6.925 


X 


10® 


2.423 


X 


10® 


2.177 


X 


10® 


BH WDXBt 


3.223 


X 


10'^ 


4.161 


X 


lO'' 


4.313 


X 


10^ 


5.088 


X 


10* 


6.099 


X 


10® 


2.843 


X 


10® 


2.920 


X 


10® 


S-Symb 


2.363 


X 


10^ 


1.946 


X 


10^ 


2.713 


X 


10^ 


2.858 


X 


lOi 


4.470 


X 


10^ 


1.378 


X 


10^ 


1.271 


X 


10^ 


D-Symb 


3.389 


X 


10^ 


3.837 


X 


10^ 


3.369 


X 


10^ 


3.369 


X 


10^ 


5.044 


X 


10^ 


2.702 


X 


10^ 


2.939 


X 


10^ 


nHe MSG 


1.497 


X 


10* 


1.017 


X 


10* 


2.116 


X 


10^ 


4.578 


X 


10^ 


4.402 


X 


10* 


1.156 


X 


10* 


1.572 


X 


10* 


gnt MSG 


4.609 


X 


10* 


7.355 


X 


10* 


4.390 


X 


10* 


1.456 


X 


10* 


1.292 


X 


10* 


3.182 


X 


10* 


5.554 


X 


10^ 


WDWD DD 


6.018 


X 


10* 


6.028 


X 


10» 


4.919 


X 


10* 


5.759 


X 


10' 


1.041 


X 


109 


1.219 


X 


108 


4.048 


X 


108 


WDNS DD 


2.725 


X 


10^ 


1.866 


X 


10^ 


4.090 


X 


10* 


1.192 


X 


10^ 


6.959 


X 


10^ 


1.800 


X 


10^ 


2.955 


X 




NSNS DD 


2.427 


X 


10^ 


2.160 


X 


10^ 


1.307 


X 


10*^ 


7.949 


X 


10* 


8.313 


X 


10® 


2.367 


X 


10® 


2.249 


X 


10® 


He DDRch 


3.076 


X 


10^ 


2.913 


X 


10' 


1.003 


X 


10' 


3.959 


X 


10^ 


1.353 


X 


10** 


2.858 


X 


10' 


2.371 


X 


10' 


GO DDRch 


2.465 


X 


10^ 


3.299 


X 


10' 


1.261 


X 


10'5 


7.225 


X 


10^ 


4.125 


X 


10' 


2.127 


X 


10' 


1.879 


X 


10' 



during the common-envelope (CE) phase that is an integral 
part of CV formation. 

As a typical example of classical CV formation in 
Model A let us consider stars of mass 2.0 and 0.2 M© in 
an orbit with initial separation 100 R© . The more massive 
star evolves to fill its Roche-lobe on the EAGB when the 
separation has fallen to TOR© (owing to tidal interaction), 
a CE forms and a HeHG star separated from the 0.2 M© 
MS star by 1.1 R© emerges. The HeHG star evolves to a 
COWD and shortly afterwards the MS star fills its Roche- 
lobe creating a classical CV. The same system evolved with 
Model B has expanded to a = 105 R© (owing to stellar wind 
mass-loss) when the first phase of RLOF begins. As a re- 
sult the post-CE binary is too wide to interact within Tgai. 
However, if the same masses are evolved from a closer initial 
separation of 50 R©, in Model A the stars coalesce during 
CE while Model B produces a classical CV. 

A typical CV Symb system can form from a 3.2 and a 
1.0 Mq star initially separated by 180 R©. RLOF first occurs 
when the primary has evolved to the EAGB and tides have 
reduced the separation to 160 R© . Mass transfer is dynami- 
cal so that a CE forms from which a close binary consisting 
of a HeHG and a MS star emerges. The helium star evolves 
to a COWD and a second phase of RLOF begins when the 



companion reaches the HG. This time the mass transfer pro- 
ceeds on a thermal timescale so that when the new primary 
reaches the GB its mass has dropped to 0.3 M© and a second 
CE phase is avoided. Stable mass transfer then continues on 
the GB until the primary loses its envelope and a DD sys- 
tem is born. This system has evolved through both a GK Per 
and a CV Symb phase. If it were to be evolved without tidal 
interaction the binary formed after the CE phase would be 
wider. The new primary would still fill its Roche-lobe on the 
HG, but only when closer to the GB, so that upon reaching 
the GB a second CE occurs in which the system coalesces. 
Thus a GK Per phase still exists but the CV Symb phase 
does not. 

Unlike the other types of CV, the sdB population in- 
creases in formation rate in Model B. As already mentioned, 
there are two main regions of the parameter space shown in 
Figure ^ from which sdB systems form. First consider the 
population formed from close orbits which can be character- 
ized by a binary of initial masses 3.2 and 1.0 M© separated 
by 20 R© . The primary fills its Roche-lobe on the HG and its 
mass has fallen to 0.6 M© by the time it reaches the GB so 
that CE evolution is avoided and stable mass transfer contin- 
ues. Eventually it loses its entire envelope to become a HeMS 
star. The companion evolves to the GB and fills its Roche- 
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Figure 5. The region of parameter space from which cataclysmic variables (CVs) form when the initial mass of one star is fixed at 
3.2 M0. Systems formed from Model A are shown as + symbols, from Model B as X symbols, with any overlap as filled squares. The 
indices i and j can assume the value of either 1 or 2 depending on which mass represents the primary star on the ZAMS, i.e. if the 
3.2 Mq star is the more massive then j = 1 and j = 2. 



lobe, by which time the HeMS star has become a COWD and 
the resulting CE evolution leaves a COWD/HeMS pair in a 
close orbit. When the new HeMS evolves to fill its Roche- 
lobe a sdB system is formed. The population formed from 
wider orbits involves a slightly more complicated path rep- 
resented by a binary with M2 = 3.2 Mq, Mi = 10.0 M© and 
a = 500 R© . A first phase of RLOF begins when the more 
massive star is on the HG with the mass transfer becoming 
dynamical when it reaches the GB. The ensuing CE leaves 
a MS/HeMS pair separated by 10 R0. After the helium star 
evolves to the helium GB it fills its Roche-lobe and mass 
transfer proceeds until all the envelope is lost and the star 
becomes an ONeWD. The MS star now begins a third phase 
of RLOF when it reaches the HG. During this phase the 
ONeWD accretes enough material to swell-up and become a 
TPAGB star so that shortly afterwards the system evolves 
into another CE from which a close HeMS/ONeWD pair 



emerges. Yet another phase of RLOF begins when the new 
HeMS star fills its Roche-lobe and the sdB system is born. 
Mass transfer continues until an AIC changes the ONeWD 
to a NS. The increase in Model B birth rate for the sdB 
population is explained by the cluster of additional systems 
about M2 = 2.6 M0 and a — 20 Rq. These evolve via a sta- 
ble RLOF phase followed by two instances of CE evolution 
which produce a close HeMS-COWD pair. When the helium 
star expands to fill its Roche-lobe the sdB phase begins. If 
tidal evolution is included the system formed after the first 
CE is much closer so that coalescence occurs during the sec- 
ond CE. 

A decrease in the common-envelope efficiency parame- 
ter (Model C, QfcE = 1) increases the CV birth rate because 
wider systems are brought within an interaction distance 
after the CE phase. Thus systems that don't interact in 
Model A can in Model C. This is countered to some degree 
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Figure 6. The region of parameter space from wliicli double-degenerate (DD) binaries form when the initial mass of one star is fixed 
at 1.6 Mq. Systems formed from Model A are shown as + symbols, from Model B as X symbols, with any overlap as filled squares. The 
unstable behaviour of some near equal-mass systems is due to the relatively short lifetime of the GB stage. In this case there is a range 
of separations where the first phase of RLOF begins with both stars on the GB and CE gives a HeWD-HeWD pair. If the mass of either 
star changes slightly then one star is still on the MS when RLOF begins and the subsequent evolution is altered. 



by closer systems that would interact in Model A now coa- 
lescing. The net effect is an increase in production. It could 
be argued that a combination of Models B and C would 
produce CV birth rates similar to those of Model A. 

Comparing the birth rates of Algols in Models A, B 
and C reveals no significant variation. This is expected be- 
cause RLOF generally begins while the primary is still on 
the MS when tides have had little effect on the orbit and 
also because a CE phase is not involved. Of the gnt MSC 
systems that form 55% of the population consist of a GB 
primary and a HG secondary while the remaining 45% are 
GB primary and GB secondary. In most cases the GB-GB 
systems have previously been GB-HG and before that were 
Algols. The general evolution path for these systems begins 
with RLOF while the primary is on the HG and the sec- 
ondary on the MS, the HG and GB lifetimes of the primary 



increase as it loses mass on the HG, mass transfer has re- 
duced qi below gcrit when the primary reaches the GB, and 
the secondary evolves to the HG and then the GB while the 
Roche- lobe filling primary is still on the GB. Of the mis- 
cellaneous RLOF systems with helium star secondaries 51% 
have a MS primary and a HeMS secondary, and 42% a HeMS 
primary and secondary. 

The birth rate of LMXBs shows an order of magnitude 
decrease when either Model B or Model C is compared to 
Model A. With tides a typical LMXB formation scenario 
begins with a primary mass somewhere in the range 10 
to 12 Mq, a low-mass companion with M2 — 1.0 Mq and 
P ~ 1 000 d. The massive primary evolves to the EAGB 
before filling its Roche-lobe and forming a CE system from 
which a fairly close HeHG-MS pair emerges. The helium star 
explodes in a supernova that leaves a NS. This may increase 
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Figure 7. The region of parameter space from which CV, Algol and X-ray binaries form when the initial mass of one star is fixed at 
3.2 M0 in Model A. 



or decrease the separation depending on the size of the ve- 
locity kick. After some time the MS star evolves to fill its 
Roche-lobe, by which time any eccentricity induced in the 
orbit by the supernova has been removed by the tides and 
a persistent NS LMXB is formed. If tides are not included 
then the orbit of the HeHG-MS pair is much wider so that 
the second RLOF phase occurs when the primary is on the 
GB, a CE forms and a DD NS-HeWD binary is produced. 
On the other hand if Oqe = 1 then the orbit of the HeHG- 
MS pair is much closer after the initial CE so that the HeHG 
star fills its Roche-lobe resulting in another CE in which the 
two stars coalesce. There is an increase in the LMXB birth 
rate for Model D, for which the secondary mass is chosen 
from the same IMF as the primary, because massive stars 
are more likely to have low-mass companions, increasing the 
number of systems that can become LMXBs. This is com- 
pensated by a corresponding decrease in the MXRB birth 
rate. All the remaining binary systems show a drop in birth 



rate for Model D because more systems form with low qi 
decreasing the likelihood of interaction. 

Interestingly the formation rate of DD systems increases 
when the tides are not used. Han (1998) finds that a phase 
of stable RLOF followed by CE evolution (RLOF-I-CE) and 
a combination of two CE phases (CE-fCE) are the two main 
channels for forming DDs with the former channel more 
likely to occur for close orbits. Adding tides to the model 
removes many systems from the closer RLOF-I-CE channel 
because the separation decreases while the secondary evolves 
on the HG and GB, after the first RLOF phase, so that the 
CE produces a much closer pair of HeMS stars. This pair 
then evolve into contact before a DD can form. Figure ^ 
shows the parameter space in the Mj-a plane for DD for- 
mation in Models A and B when the initial mass of one 
star is fixed at 1.6 Mq. The systems with M\ > 8.0 Mq are 
WD-NS binaries that survive the supernova kick while the 
range of close systems with a <, 100 Rq are RLOF-I-CE WD- 
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Figure 8. The region of parameter space from which CV, Algol and X-ray binaries form when the initial mass of one star is fixed at 
3.2 M0 in Model F. 



WD binaries, of which there are noticeably more formed in 
Model B. A typical RLOF+CE WD WD DD formation sce- 
nario in Model B begins with Mi = 2.2 Mq, A'h = 1.6 Mq 
and P = 2d. The primary fills its Roche-lobe on the HG 
when P = 1.9 d. RLOF ends on the GB just before the 
primary becomes a CHeB star, loses what little is left of 
its envelope, and moves to the HeMS. Now Mi = 0.35 Mq, 
M2 = 3.45 M© and P = 60d. The orbit has widened as a 
result of the conservative mass transfer. The more massive 
star (now assuming the role of the primary) evolves to the 
GB and overflows its Roche-lobe on a dynamical timescale 
so that a CE forms from which a HeMS-HeMS pair emerges 
with P = 0.1 d. Both of these evolve to WDs in a close 
COWD-COWD DD system. If the same initial parameters 
are evolved with tides the system has P = 30d just before 
the CE so that a closer HeMS-HeMS pair forms. This evolves 
into contact before either star can become a WD. If the same 
masses are evolved with a longer initial period, P ~ 240 d, 



then a DD forms via the CE-I-CE channel in Model B. In 
this case the primary fills its Roche-lobe on the EAGB, CE 
evolution ensues and a HeHG star with AIi — 0.6 Mq in an 
orbit of P = 6d with the MS star emerges. The HeHG star 
evolves to a COWD. The MS star later fills its Roche-lobe 
on the GB, leading to another CE and the formation of a 
HeWD-COWD DD system with P = 0.03 d. Evolving this 
case with Model A leaves a closer binary after the first CE 
so that the second CE ends in coalescence of the two stars. 
Systems with wider initial orbits can survive the second CE 
phase to produce a DD binary when either Model A or B 
is used. Even wider systems form a WDWD DD with no 
interaction but these have much longer periods than those 
formed via CE and are not easily detected. 

Inspection of Figure ^ reveals a large empty region in 
parameter space between the close group of systems that 
form DDs via RLOF-I-CE and the wider group of systems 
that form DDs via CE-I-CE. In this gap the initial period 
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Figure 9. The probability, as a function of time, of a binary being found as a member of certain individual populations of Model A. 
We assume that all stars are born at the same time. The individual populations shown are: blue stragglers (full line), cataclysmic 
variables (dashed line), Algols (dash-dot line), persistent X-ray binaries (dotted line) and double-degenerates (dash-dot-dot-dot line). 
Each population includes contributions from each of its sub-populations listed in Table ^. The numbers for double-degenerates have been 
reduced by a factor of 10 and the numbers for persistent X-ray binaries have been increased by a factor of 10. 



is short enough that the second CE leads to coalescence in 
both Models A and B. However pockets of DDs do exist in 
Model A because stable RLOF can occur on the HG. Mass 
transfer reduces the primary mass and actually causes a net 
increase in the separation. This phase of RLOF, coming after 
the first CE, is non-conservative so that under the influence 
of tidal synchronization the separation initially decreases. 
As the primary mass is reduced further and the strength 
of the tides weakens the orbit starts to expand. When the 
primary reaches the GB a second CE occurs but the orbit is 
wide enough, and the size of the primary small enough, that 
coalescence is avoided and a DD formed. This is however a 
very unstable region of phase space. If the initial separation 
is reduced slightly the RLOF phase after the first CE oc- 
curs on the MS, the primary mass falls until it develops a 
convective envelope when dynamical mass transfer leads to 



coalescence. A slightly larger initial separation means that 
the RLOF phase starts on the GB before the star has lost 
any mass, so it is much larger and the CE ends in coales- 
cence. If Ml is increased slightly the RLOF phase after the 
first CE still starts on the HG but, owing to a larger gi, leads 
to dynamical mass transfer and coalescence. A reduction in 
Ml means the primary takes longer to fill its Roche-lobe be- 
fore the first CE so that the tides give a greater reduction in 
the separation, the system formed after the CE is closer, and 
the next RLOF phase starts on the MS, once more ending 
in coalescence. With Model B these regions do not occur at 
all because the orbit is always closer when the second CE 
forms so that the binary never survives. 

The lower metallicity population. Model E, has in- 
creased birth rates in almost all cases compared with 
Model A. This behaviour can be primarily attributed to 
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the shorter nuclear lifetimes of lower metallicity stars with 
M ^ 9 Mq (see Figure 5 of Papl). For example, a solar mass 
star with Z = 0.0001 evolves to a WD cooling cooling track 
in half the time it takes a solar mass star with Z = 0.02. 
This gives it more chance of interacting with its compan- 
ion when in a relatively wide orbit. Model E represents a 
very young population that would have formed soon after a 
galaxy condensed. 

It should be noted that the choice of separation distri- 
bution ^ (In a) and its limits, particularly the upper cut-off, 
affects the calculated birth rates and numbers. Only the BSS 
and D-Symb populations generate any members from initial 
separations greater than lO** Rq but these are rare and so 
do not affect the absolute numbers very much. 

4-2.3 The Role of Semi-Latera Recta 

Figures ^ and ^ show the parameter space of CV, Algol and 
X-ray binaries that are formed for Models A and F respec- 
tively when the initial mass of one star is fixed at 3.2 M©. 
In Model F we form binaries with eccentric orbits and, as 
expected, this means that systems with wider orbits can in- 
teract at periastron and contribute to the rates. This does 
not automatically lead to an increase in formation rates be- 
cause closer systems are more likely to be destroyed by a 
collision or coalescence. If we plot semi-latus rectum on the 
vertical axis of Figure ^ then it appears virtually identical 
to Figure |^ Therefore, if a distribution of orbital angular 
momentum or I is used to determine the initial state of each 
binary, rather than one of semi-major axis or period, the 
results do not depend on the form of any chosen eccentric- 
ity distribution. In fact eccentricity need not be a free pa- 
rameter. This becomes evident when considering that tidal 
interaction conserves angular momentum and that almost 
all systems circularize before RLOF. Therefore (see eq. ^ 
systems with the same initial / end up in a circular orbit 
with the same separation when RLOF starts, and their sub- 
sequent evolution will be identical. The only instances where 
the initial eccentricity distribution does matter are binary 
populations in which tidal interaction is weak, such as wide 
binaries that may avoid RLOF (e.g. the Barium stars, see 
Karakas, Tout & Lattanzio 2000) 

The circularization method of Model G usually has in- 
creased formation rates because it acts to bring systems 
closer making them more likely to interact. Unlike Model F 
none of the angular momentum is absorbed by the stars so 
that systems do not get so close and as a result less are 
destroyed. 

The relative fractions of blue stragglers, cataclysmic 
variables, Algols, persistent low-mass X-ray binaries and 
double-degenerates as a function of time for Model A are 
shown in Figure ^, assuming that all binaries were born in 
a single burst of star formation. This represents the proba- 
bility of a binary observed in a star cluster being of a par- 
ticular type. Each population is a combination of its sub- 
populations as listed in Table |^ The peak in the blue strag- 
gler distribution at approximately 5 000 Myr is primarily due 
to a decrease in progenitor systems and its timing depends 
on the ratio of the mean binary separation, which remains 
constant, to the radius of the largest MS star, which de- 
creases in time (Hurley et al. 2001). Both the CV and DD 
probabilities show a steady increase in time, reflecting both 



the long-lived nature typical of these systems and the con- 
stant production of WDs from the evolving mass distribu- 
tion. As the existence of LMXBs is directly linked to the 
production of NSs or BHs it is no surprise that their prob- 
ablity distribution peaks at an early age. The probability 
of observing an Algol also shows a peak that is linked to 
the evolution of massive stars as these are more likely to 
fill their Roche-lobes on the MS or HG, thus avoiding a 
common-envelope phase. 

4.3 Comparison with Observations 

To determine whether the parameters of Model A are a good 
choice we can compare with observations. For most of the 
individual binary populations observational birth rates or 
numbers in the Galaxy are very uncertain because of se- 
lection effects involved when undertaking surveys. However 
enough data exist overall to enable a meaningful compari- 
son with the results from which a decision on the best model 
parameters can be made. We should add that caution must 
be exercised when comparing observationally defined classes 
of binary with the model classes defined in this work. Some 
overlap between classes necessarily exists: owing in part to 
the restrictive nature of the theoretical model and to remain- 
ing uncertainty about the details of some of these objects in 
the first place. 

4-3.1 Cataclysmic Variables 

Ritter & Burkert (1986) estimate the CV birthrate in our 
Galaxy to be about 10"" pc"^ yr"^ or 6 x 10"^ yr"^ if 
Vgai = 6 X lO^^pc^. In calculating this rate they corrected 
the observed local space density for the selection effect of 
only a small fraction of CVs being bright enough to be seen 
even in a distance limited sample. They assumed that the 
mass spectrum of WDs in CVs is not too different from that 
of single WDs. The birth rate of classical CVs in Model A 
is 2.0 X 10~^ yr~^ which is only a factor of 3 greater than 
this so that the two are consistent, considering that a bi- 
nary fraction less than unity would reduce the derived rate. 
The intrinsic local space density of CVs quoted by Ritter & 
Burkert (1986) is 1 - 2 x 10"" pc"^ (6 to 12 x lO'' in the 
Galaxy), which they note seems rather high, and can only 
be matched by the results of Model D. Yungelson, Livio & 
Tutukov (1997) find the present local space density of clas- 
sical CVs in their model to be about 3 x 10~^ pc~'^ which 
is in agreement with the Hertz et al. (1990) observed num- 
ber and consistent with the results of all the models except 
Model D. 

4.3.2 Algols 

Observational data on Algol systems suffers from the lack 
of a comprehensive survey aimed at establishing the local 
space density. The catalogue of eclipsing binaries published 
by Brancewicz & Dworak (1980) contains parallaxes for 653 
systems classified as Algol. Apparent magnitudes for these 
systems can be obtained from the catalogue of interacting 
binaries published by Wood et al. (1980) although no men- 
tion is made as to the completeness limit of this catalogue. 
Figure ^ shows the distribution of the Algol systems that 
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classical Algols with Mi > M© currently in the Galaxy is 
then 1.2 x 10^. However, almost half of the systems classi- 
fied as Algol by Brancewicz & Dworak (1980) are assessed 
to have an extremely unlikely probability of RLOF by Bud- 
ding (1984) who analysed the semi-detached nature of each 
system. Therefore it is likely that this number is an over- 
estimate. On the other hand, only about 1/3 of all classical 
Algols are likely to be observed as eclipsing binaries (see 
Appendix A2 for details). The number of Algols with a MS 
or sub-giant primary more massive than 1 Mq is predicted 
by Model A to be about 1.4 x 10^ (combining the MS, cold 
and hot Algol populations). This is an order of magnitude 
greater than the number predicted by observations but is 
close enough to not be a major worry when we consider the 
uncertainties. 



Figure 10. The apparent magnitude distribution of Algol sys- 
tems identified by Brancewicz & Dworak (1980). The magnitudes 
are taken from Wood et al. (1980). A best fit to tlie data for 
m < 9.5 is shown by the full line. 



have a primary mass greater than 1 Mq as a function of their 
apparent magnitude. The apparent brightness of a star is 
linked to its distance and bolometric luminosity by 

^ " 4^ 

which defines a volume 



(107) 



V = 



47r (J^Y^^ 
3 UttJ 



(108) 



Assuming that the stars are uniformly distributed within 
this volume the number of stars out to apparent brightness 
I is 



Ni (m) = Ar^''^n(m) 



(109) 



where A is some constant. The function n(m) depends only 
on the stellar mass and incorporates such quantities as the 
mass function, the dependence on luminosity and how long 
a star of a certain mass remain s vis ible. Integrated over all 
masses this is constant in eq. ( [l09| ) so it acts only to nor- 
malize the relation. 

To determine the completeness limit of the Alg ol sam- 
ple we made a series of least-squares fits of eq. (109) to the 
data starting from the lowest magnitude bin, each succes- 
sive fit includes data from the next bin in line. The number 
of Algol systems in the Brancewicz & Dworak (1980) sam- 
ple that are brighter than 10th magnitude is 282. A least- 
squares fit to this subset of the data is shown in Figure |l^ 
and gives 300 systems brighter than 10th magnitude which 
is just consistent with the data within the Poisson error. 
If we include systems above 10th magnitude in the fit this 
leads to a large increase in the RMS error between the fitted 
function and the data as well as making the two inconsis- 
tent within Poisson error. Therefore it seems safe to assume 
that the Algol sample is complete to 10th magnitude. This 
corresponds to a solar mass primary at the end of the MS 
being visible to a distance of roughly 180 pc. Returning to 
the Brancewicz & Dworak (1980) sample and counting only 
systems within this distance with primary mass greater than 
1 Mq we find 50 Algol systems corresponding to a local 
number density of n ~ 2 x 10~^ pc~^ and the number of 



4.3.3 Low-Mass X-ray Binaries 

The LMXB catalogue compiled by van Paradijs (1995) lists 
about 80 bright persistent sources in the Galaxy, making 
no distinction as to whether the compact star is a NS or 
BH. This number is in good agreement with the results of 
Model A, and in fair agreement with Model B. The birth 
rate of persistent NS LMXBs in Model A is similar to the 
rate found by Portegies Zwart & Verbunt (1996) in their 
standard model that includes velocity kicks. The relative 
distribution of NSs and BHs in X-ray binaries depends on 
the choice of single star mass abov e which BHs are formed 
rather than NSs (see Section 2.6.6). This mass is not well 
constrained. It also depends on the assumption that a NS 
collapses to a BH when it accretes enough material to take 
its mass above 1.8 Mq (Bombaci 1996; see also Papl). If 
the Eddington limit for mass transfer is applied then the 
birth rate of persistent NS LMXBs in Model A increases to 
3.612 X 10"® while the BH LMXB rate drops to 1.446 x 10"^ 
The actual number of persistent LMXBs drops in this case, 
making the model inconsistent with the observations. 



4.3.4 Symbiotic Stars 

About 125 symbiotic stars are observed, with WD sec- 
ondaries about 10 times more prevalent than MS secon- 
daries (Kenyon 1986). Estimates of the total number in 
the Galaxy range from about 1 000 (Boyarchuk 1975) up 
to 20 000 (Luthardt 1992). Of the known symbiotic stars 
roughly 20% are D-Symb but this is probably biased be- 
cause these are harder to identify due to dust enshrouding 
the system. Kenyon (1986) predicts from theory that the 
total number in the Galaxy should be about 3 000 with S- 
types and D-types contributing equally, although the rela- 
tive number depends on the period distribution. Yungelson 
et al. (1995) consider only symbiotics with WD secondaries 
and derive a birth rate of 0.073yr~^ with 3 370 present in the 
Galaxy using their standard model. These numbers are in 
good agreement with the predicted numbers from all of the 
models except Model D which has only about 350. At least 
98% of the symbiotics found in each of the models have WD 
secondaries, the remainder have HeMS secondaries, and the 
birth rates don't appear to fluctuate significantly between 
models, except Model D. 

As already mentioned, the systems with MS secondaries 
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are part of the cold Algol population, the number of which 
is typically 7 x 10^ when only GB or AGB primaries are 
considered. This is much larger than the expected number so 
perhaps they are only observed as symbiotic for some of the 
RLOF phase or possibly only a small fraction are observed 
as symbiotic at all. The birth rate is typically of the order 
10~^yr~^ so an average symbiotic lifetime of 10* yr for these 
MS+giant systems would be consistent with the predicted 
numbers. The large CV Symb numbers are not consistent 
with these systems being observed as symbiotic. It could 
be that they are simply not found, because nova explosions 
don't occur on the WD surface, and the giant primary most 
likely dominates any light from the accretion disk. 

4-3.5 Double-Degenerates 

Maxted & Marsh (1999) have conducted a radial velocity 
survey to measure the fraction of DDs among hydrogen-line 
(DA) WDs. They claim that their survey is the most sensi- 
tive so far undertaken, in terms of the detection of DDs, yet 
they only found two systems in a sample of 46 DA WDs, and 
these two were known to be DD prior to the observations. 
Their detection efficiency is 80% or higher for periods less 
than 10 d and quickly drops to zero for wider orbits. Thus 
it seems that the fraction of DDs with P < 10 d among DA 
WDs has a lower limit of i/wd ~ 0.04. This is consistent 
with the survey of Bragaglia et al. (1991) who find between 
one and three DDs in a sample of 49 DA WDs. Combining 
their results with a theoretical model that predicts the pe- 
riod, mass and mass-ratio distributions of DDs, Maxted & 
Marsh (1999) find a 95% probability that the DD fraction 
lies in the range 0.017 to 0.19, independent of the details of 
the model used. Saffer, Livio & Yungelson (1998) conducted 
a radial velocity survey of a sample of 107 DA WDs and 46 
sdB stars. The sdB, or helium, stars were included because 
they are believed to evolve to the WD cooling track on a rela- 
tively short time, following the exhaustion of central helium. 
From this sample they found a possible 23 binaries of which 
14 are WD-HWD, 7 sdB-hWD and 2 WD-hMS. Their detec- 
tion efficiency is only 70% for periods between 2 and 20 d. 
Of the 14 possible WD-I-WD binaries, seven have now been 
confirmed as DDs (Marsh, private communication) from the 
sample of 107 DA WDs so this suggests Uvju > 0.07. 

Consideration of all these results together requires a 
reasonable theoretical model to predict at least Uvju ~ 0.05 
for DDs with periods less than 10 d, and z^wd « 0.1 when all 
DDs are taken into account. The frequency of DDs among 
DA WDs can be converted to a DD birth rate if it is as- 
sumed that one DA WD is born per year (Fleming, Liebert 
& Green 1986; Phillips 1989): equivalent to our one star of 
M > 0.8 M© per year. In terms of the WDWD DD birth 
rates both Models A and B satisfy the requirements of the 
observations: they predict 0.113 and 0.123 yr~^ respectively 
when all periods are included; if only DDs with P < 10 d are 
considered the rates drop to 0.053 and 0.059 yr"'^. Model C 
is in substantial disagreement with the observations: it pre- 
dicts 0.076 yr~^ for all periods and only 0.013 yr~^ for 
P < 10 d; Model D is even worse. The predicted rate when 
eccentric orbits are allowed is also low but the small period 
rate in this case is still 0.051 yr~^, in good enough agree- 
ment with both Model A and the observations. Han (1998) 
includes periods up to at least 100 d in his calculation of 
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Figure 11. The distribution of periods of WDWD DD systems at 
formation in Model A (full line). The same distribution weighted 
by the gravitational radiation timescale of each system is shown 
for Model A (dashed line) and Model B (dotted line). Also shown 
are the periods of 22 observed DD systems (asterisks) provided by 
Marsh (private communication) Note that the vertical positioning 
of these points bears no relation to any physical parameters. 



0.03 yr^^ for the DD birth rate but his model uses Oce ~ 1 
(although Han accounts for the ionization energy in calcu- 
lating the envelope binding energy so this corresponds to a 
higher and variable Qqe in this work). 

The period distribution at formation for DDs in 
Model A is shown in Figure ^ It has a peak at about 0.05 d. 
Included are periods for 22 suspected DD binaries provided 
by Marsh (private communication) . To account for the short 
lives of short-period systems we also show the period distri- 
bution weighted by the corresponding gravitational radia- 
tion timescale for each system. The weighted distribution 
for Model B shows that the observed distribution is much 
better represented when tidal evolution is included in the 
model. 

If the merging of two WDs is the progenitor of all type 
la SNe then it is necessary for an observed population of DDs 
to contain a certain proportion of very close systems with the 
combined mass of the components near or above the Chan- 
drasekhar limit. The absence of these systems has proven a 
major obstacle for the proponents of the WD merger pro- 
genitor model (e.g. Saffer, Livio & Yungelson 1998). How- 
ever, the merger timescale decreases with shorter periods 
and larger binary masses so, as Figure ^ shows, such sys- 
tems are less likely to be observed. The statistics of DD 
searches must be improved before the WD merger progenitor 
model can be ruled out on this basis. In Model A the percent- 
age of DDs with P < 10 d present in the Galaxy now that 
contain two COWDs, have Mi -I- M2 > Mch and will merge 
within the lifetime of the Galaxy is 0.3. The percentage con- 
taining a COWD and a HeWD of mass at least 0.15 M© that 
will merge is 2.3. If all periods are included then the per- 
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centages drop by a factor of 10. Considering that 30 DD 
systems with measured periods are known (Marsh, private 
communication) we predict that a possible type la SNe pro- 
genitor should soon be discovered (e.g. Maxted, Marsh & 
North 2000). 

The total number of DDs currently in the Galaxy is 
predicted by Model A to be 6x 10*. If WDs are detectable for 
only the first 10* yr after formation then the number drops 
to 10^, more than the 3 x 10^ predicted by the standard 
model of Han (1998). 

Bragaglia, Renzini & Bergeron (1995) found that, of 
a sample of 164 DA WDs, 15 were low-mass WDs with 
M < 0.4 M0. Low-mass WDs are interesting because they 
are not predicted by standard theory of single star evolution 
so they are believed to be, or have been, members of close bi- 
naries. They are not however directly representative of the 
number of DDs because they can be produced via other 
channels of binary interaction, such as in WD/NS binaries. 
Our theoretical models include WDs with A/ < 0.5 Mq in 
the LMWD populations, so a birthrate of 0.1 yr^^ is a lower 
limit for comparison and is in good agreement with all mod- 
els except Model C. 



4.3.7 r -ray Bursts 

The observed rate of 7-ray bursts is about lO"*" yr^^ per 
galaxy (Piran 1996) which is typically an order of magnitude 
less than the merger rate of NS binaries found in all models 
except Model D. Numerical models of NS mergers suggest 
that the 7-rays produced are beamed and appear only in 
certain directions (Davies et al. 1994). An optically thick 
radiation-electron-positron plasma, or fireball, with initial 
energy larger than its rest mass is formed. The fireball is 
highly aspherical and expands along the polar axis to form a 
jet. Only when the fireball reaches ultra-relativistic velocities 
and its material becomes optically thin, or reacts with the 
interstellar medium, do the observed 7-rays emerge. If the 
width of the jet is Q then we observe 7-ray bursts only from 
a fraction 29"^ of NSNS mergers (Piran 1996). Portegies 
Zwart & Yungelson (1998) find the neutron star merger rate 
in their model to be about 2 x 10~^ yr~^ which they claim is 
consistent with the observed rate of 7-ray bursts if the 7-ray 
emission is beamed into an opening angle of 10°. The same 
reasoning can reconcile the rates for all our models, except 
Model D, with the observed rate of 7-ray burst events. 



4.3.6 Supernovae 

Cappellaro et al. (1997) combined the results of five inde- 
pendent SN searches to obtain a sample of 110 SNe from 
which they derived the following birth rates for the Galaxy: 

4±lxl0"^ yr"^ SNla 
2 ± 1 X 10"^ yr"^ SNlb/c 
12 ± 6 X 10~^ yr-i SNll. 

The SNII rate is consistent with all the Z = 0.02 models 
except possibly Models C and G. Note that the increase in 
type II SNe for the lower metallicity population (Model E) is 
due to supernovae occurring in lower mass single stars (see 
Table 1 in Pols et al. 1998). An obvious discrepancy lies in 
the birth rates of SNIb/c because all the theoretical rates are 
too high, except for Model D which is in good agreement. 
This can be reconciled by assuming a binary fraction smaller 
than unity. 

The SNll rate is a fairly robust rmmber but the SNIb/c 
rate is sensitive to many of the assumptions underlying the 
model. In all models other than Model D the majority of 
SNIb/c come from naked helium stars formed from relatively 
low-mass progenitors (10 — 2OM0) stripped by binary inter- 
action. In a model containing only single stars with Z = 0.02 
the SNII rate is the same as for Model A but the SNlb/c rate 
is reduced to 2.7 x 10~^ yr~^ which would be consistent with 
observations. These single SNe Ib/c come from stars ini- 
tially more massive than 24 M0 which become Wolf-Rayet 
stars prior to explosion. However this rate is sensitive to 
the assumed mass-loss rate for red supergiants (see Papl), 
as well as to our assumption that BH formation ignites a 
supernova explosion, both of which are very uncertain. If 
mass-loss rates are lower or if black holes form without a 
SN, the SNlb/c rate from single stars would be lower. 



5 CONCLUSION 

The primary purpose of this paper was to describe in detail 
a rapid evolution algorithm for binary stars. Recognizing 
that tidal interaction is an important process in binary evo- 
lution we have also studied its effect on the results of binary 
population synthesis and made a global comparison with 
observations. 

The main effects of tides on the evolution of an inter- 
acting binary are: 

(a) circularization of the orbit, generally well before RLOF 
occurs, and 

(b) the exchange of angular momentum between the orbit 
and the spins of the components. 

Because as a star expands tidal synchronization transfers or- 
bital angular momentum to the stellar rotation, in general 
the latter effect tends to bring binary components closer to- 
gether and cause them to follow an evolutionary path similar 
to that of a closer binary if tides were ignored. However, we 
find that as a result of the sometimes very convoluted evolu- 
tionary paths of interacting binaries, this general principle 
does not always hold and it is difficult to summarize the 
overall effect of tides on binary evolution. 

The primary differences in our population synthesis re- 
sults when tidal evolution within binary systems is ignored 
are as follows: 

(i) a 35% decrease in the birth rate of classical CVs but 
a 23% increase in sdB CV production; 

(ii) lower birth rates for LMXBs with the persistent NS 
LMXB rate decreasing significantly but conversely the num- 
ber of NS LMXB increases; 

(iii) a 10% increase in the birth rate of double- 
degenerates; 

(iv) an increase of 190% in the incidence of exploding 
HeWDs; 

(v) and 55% more Algols currently in the Galaxy. 
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While these difFerences may seem substantial they are gen- 
erally not enough to confirm one way or another whether 
tidal evolution helps to explain the observed binary pop- 
ulations. This is mainly due to the lack of comprehen- 
sive binary searches and the selection effects involved that 
make discrimination on the basis of observational evidence a 
risky business. However, tides are present in binary systems 
whether we like it or not. One area where neglecting tidal 
evolution does seem in clear conflict with observations is in 
the incidence of HeWDs that explode as supernovae. This 
rate is too high by a factor of two to be accounted for by 
either the observed numbers of type la or type lb SNe (see 
Tout et al. 2001 for further discussion). 

The common-envelope efficiency parameter Qce is an 
uncertain factor in a phase of evolution that is crucial for 
the production of many types of binary. The failure of the 
model with Oqe = 1 to produce anywhere near enough DDs 
would seem to favour acE = 3. Choosing the secondary mass 
independently of the primary mass does not produce enough 
DDs, symbiotics or type la SN candidates. All in all the 
results favour using the properties of Model A. In fact the 
standard tidal model, represented by Model A, is not in 
disagreement with any of the observations, except in the 
production of too many type Ib/c SNe. This discrepancy 
could indicate a lower mass loss rate for helium stars. 

A major conclusion of this work must be that it is ex- 
tremely difficult to set contraints on any of the parameters 
involved in binary evolution from population synthesis of 
birth rates and galactic numbers of the various types of bi- 
nary. This is certainly true while such a large number of 
parameters remain uncertain. The task would be simpler 
if we could find observational tests specific to an individ- 
ual parameter. For example, Algols do not require common- 
envelope evolution for formation. Therefore, if observational 
constraints on the number of Algol systems currently in 
the Galaxy improve then this could provide a suitable test 
for models of thermal-timescale mass transfer. The effect 
of varying additional model parameters, such as the binary 
enhanced mass loss, and initial conditions, such as the sepa- 
ration distribution, would need to be quantified before such 
a test could be reliably implemented. This is beyond the 
scope of this paper but will be the subject of future work. 
In terms of tides it is difficult to isolate direct tests of the 
tidal evolution model with the type of population synthe- 
sis performed in this paper. In fact, to properly constrain 
the strength of tidal interaction it is necessary to study pre- 
RLOF binaries, for example, as has been done previously by 
Zahn (1975, 1977). What we would like to emphasise is that 
the outcome of binary evolution is sensitive to the phys- 
ical processes of tidal circularization and synchronization. 
Therefore, any attempt to constrain uncertain parameters 
in binary evolution by the method of population synthesis 
must utilise a binary evolution algorithm that incorporates 
a working model of tidal evolution, such as that presented 
here. 

An exciting challenge for the future involves attempt- 
ing to reproduce the individual orbital characteristics of a 
large number of observed binaries with our binary evolution 
model. This will allow multi-parameter fitting and should 
become a powerful tool as the statistics of binary surveys 
continue to improve. Our preliminary work in this area has 
shown that the observed properties of DD binaries are much 



easier to explain when tidal evolution is included. A detailed 
exploration of the DD parameter space will be the focus of 
another paper. 

Tidal synchronism is important but because orbits gen- 
erally circularize before Roche-lobe overflow the outcome of 
the interactions of systems with the same semi-latus rectum 
is almost independent of eccentricity. Although the inclusion 
of a distribution of eccentricities seems natural it is not nec- 
essary in population synthesis of interacting binaries, how- 
ever, the initial separations should be distributed according 
to the observed distribution of semi-latera recta rather than 
periods or semi-major axes. 

A necessity for the near future is a thorough compar- 
ison of our BSE algorithm with the workings of a detailed 
evolution code (e.g. Nelson & Eggleton 2001). In this way 
we can improve the algorithm, especially the treatment of 
the Hertzsprung gap phase, and hopefully add more strin- 
gent constraints to many of the evolution variables. Work 
is already underway to provide mo re acc urate descriptions 
of t he pa rameters k'2 (see Section 2.3.1) and A (see Sec- 



tion 



2.7.1), through investigation of the detailed stellar mod- 
els provided by Pols et al. (1998). The algorithm will also 
be improved by providing options for how angular momen- 
tum is lost from the binary system during non-conservative 
mass-transfer (see Section 2.6), reflecting the various possi- 



ble modes of mass-transfer (Soberman, Phinney & van den 
Heuvel 1997). When this work is completed we will per- 
form additional population synthesis calculations in order 
to quantify how the various improvements affect the results 
presented here. 



6 AVAILABILITY OF THE BSE CODE 

The BSE algorithm allows the entire evolution of even the 
most complicated binary systems to be modelled in less than 
a second of CPU time rather than the several hours required 
when using a full stellar evolution code. It is therefore ideal 
for synthesising large populations of binary stars and will 
prove an extremely useful tool for testing uncertain pro- 
cesses in binary evolution. To obtain a copy of the BSE 
package described in this paper send a request to the au- 
thors who will consider providing the FORTRAN subroutines 
by ftp. 
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APPENDIX A: 

Al Orbital Parameters After Supernova 

Consider a frame of reference in which the pre-supernova 
centre-of-mass is at rest, the primary star is about to ex- 
plode and the secondary is at the origin. The initial orbit 
is in the XY-plane, so that the initial specific angular mo- 
mentum vector is directed along the positive Z-axis, and the 
separation vector r = ri — r2 is directed along the positive 
Y-axis, as shown in Figure The relative velocity of the 
stars is 

V = -f orb (sin f3i + cos (Al) 
where /3 is the angle between r and v and is such that 

1/2 



sin P 



cos P 



r (2a — r) 
e sin E 



(1 - e2cos2E)^/^ 
and E is the eccentric anomaly in Kepler's equation 
M = E — e sin E 



(A2) 
(A3) 

(A4) 



for mean anomaly A4 , which varies uniformly with time be- 
tween and 27r. The orbital speed is defined by 



vii,^r^ + r^^ = GMb(--- 
\r a 



(A5) 



As well as losing an amount of mass AMi the primary is 
subject to a kick velocity during the supernova explosion 
so that 



Ml - 
Mb - 
vi - 
where 



M[ ■ 

Ml, 
f 

V 1 



■- Ml - AMi 
: Mb - AMi 

Vi + Vfe 



Vfe = Vk (cos Lj cos (jii + sin ui cos -f sin (/>k) 



Vk = Ivfcl .(A6) 



Here i, j and k are unit vectors in the X, Y and Z directions 
respectively. We assume the separation is constant as the 
explosion is instantaneous. To find the separation at the 
moment of explosion we randomly choose a mean anomaly 
M and then solve eq. (A4) for the eccentric anomaly E using 
a Newton-Raphson method. Then 



r = a (1 — e cos E) 



(A7) 



in terms of the initial semi-major axis and eccentricity. This 
is necessary because the binary is evolved using values aver- 
aged over an orbital period so that the exact separation at 
any one time is not known for an eccentric orbit. 
After the supernova the new relative velocity is 



V„ = V + Vfe 

= (ufe cos LO cos ( 

(ufe sinoj cos q 



- Vovh sin (3) i -f 
Worb cos /3) j -f Ufe sin 0k . 



From eq. (A5) it must be true that 

2 



\r an J 



(A8) 



for the new orbital parameters, where 
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Figure Al. The binary and supernova Icick geometry 



= t^fc + t'orb ~ ZWorbffc (cOS LO COS (j) Sin /3+ 



(A9) 



sin uj cos (j) cos /3) . 

This can be solved for the semi-major axis a„ of the new 
orbit. The specific angular momentum of the new system is 

h' = r X v„ 



(AlO) 
(All) 



so it follows from eq. (^) that 
GA/;a„ (l - ef,) =lr X v„|' 
where 

|rx v„|^ = \j}f.svc?(l) + {vk cos ci; cos — Voih sin .(A12) 

This we solve for the new eccentricity e„ of the orbit. If 
either On < or Cn > 1 then the binary does not survive 
the kick. 

The angle v between the new and old angular momen- 
tum vectors is given by 

f orb sin fi — Vk cos UJ cos (f) 



^2^ 1/2 



(A13) 



[f ^sin^(^ -I- {vk cos Lu cos — Worb sin 

An amount of mass AMi is ejected from the primary, and 
hence from the system, so that the new centre-of-mass has 
a velocity 

AM1M2 



relative to the old centre-of-mass frame. This agrees with 
Brandt & Podsiadlowski (1995) for initially circular orbits. 
Note that the (f) used here is tt — as defined in their coor- 
dinate system. 

To determine the kick velocity it is necessary to choose 
Vk, the magnitude of v^, 0, the angle between and the 
orbital plane, and u, the angle between the projection of 
Vfc on to the orbital plane and the X-axis, from appropri- 
ate distribution functions. We take the kick speed from a 
Maxwellian distribution 



2 vi 



(A15) 



with a dispersion of ak ~ 190kms~^ (Hansen & Phinney 
1997) based on analysis of various pulsar proper motion sam- 
ples. 

For a variable X uniformly distributed between and 
1 we may write 



P (life) dvk = dX 

X= P(vk)dvk =F{vk) 



The change of variable 



gives us 

2 

F (u) = Erf (u) j=ue 



(A16) 
(A17) 

(A18) 
(A19) 



so that X = F (u) can be solved for some X, by Newton- 
Rapshon iteration, to give Vk- 

The kick direction is uniform over all solid angles so 

that 

P(0)=cos0 , ^<<^<|, (A20) 
and UJ is uniformly distributed between and 2tt. 



A2 Eclipsing Semi-Detached Binaries 

Consider a semi-detached system in which a primary of ra- 
dius R is filling its Roche-lobe and is separated from a com- 
panion of radius r by a distance a. The angle that a line 
which is the inner tangent to the surface of both stars makes 
with the line joining the centres of the stars can be given by 

R + r 

sm',p=^-^. (A21) 

If we assume that half of the companion must be eclipsed for 
the light-curve to be noticeably changed and that R ~ Rl, 
then 



(A22) 



Define i as the inclination angle between the line-of-sight 
and the normal to the orbital plane of the binary. If i = 
then the orbital plane is perpendicular to the line-of-sight 
and the binary will never eclipse; if j = 7r/2 then it is edge- 
on and will always eclipse. The probability of i between i 
and i + di is 



f{i) di = k sin i di , 



(A23) 



for a normalization constant k. Only one orientation of the 
orbit gives i = so it is the least likely case, whereas a 
full circle of orientations gives i — -k 12. The normalization 
constraint 

»7r/2 

/(i) di = 1 (A24) 



determines k = \. Now the binary eclipses if 

7r/2 - 7/; < i < 7r/2 (A25) 

where 

7r 

7/; — arcsin(l/3) = — — arccos(l/3) . (A26) 

Therefore the fraction of systems which are likely to be 
eclipsing is 

P = / smidi = 1/3. (A27) 

J arccos(l / 3) 



© 2000 RAS, MNRAS 000, 



